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Abstract 


Contrary to older descriptions, the air currents in the middle and upper troposphere over 
the low latitudes are not steady but highly variable. There are trains of well-developed 
waves or vortices. When these disturbances combine with troughs in the polar westerlies, 
the resulting extended troughs often reach from pole to equator. Poleward flow of heat 
must take place mainly along these troughs. Energy from low latitudes is injected into 
the temperate zone in narrow and variable areas only. 

A description of the tropical atmosphere in the layers below and above the bases of the 
cumuli leads to an attempt to explain the steadiness of the low-level wind field which 
contrasts with the high variability of the upper layers. Then the strength of the mean 
meridional circulation and the fields of divergence and convergence at the surface are cal- 
culated, followed by a summary of available information on the high-level circulation. The 
field of motion aloft varies greatly with longitude, so that inferences concerning the tropics 
as a whole can only be made from maps covering most longitudes. Such maps (300-mb charts) 
are available for a limited period down to latitude 10° N. They indicate that regions where 
the easterlies decrease upward alternate with regions of increasing easterlies. 


Introduction 


For many years, the belief has held sway steady as those near the ground. Therefore, 
that the presence of atmospheric disturbances with the aid of few and scattered upper-air 
is peculiar to middle latitudes, whereas steady data it was held possible to deduce the correct 
conditions prevail in the arctic and in the circulation picture aloft. On any given day 
tropics. This conviction helped to reduce the the upper flow should be indicative of average 
excessive number of variables which confront conditions, in contrast to middle latitudes. 
the meteorologist at every turn. In the tropics, . Solar energy accumulated in the tropics ulti- 
it found strong support in the apparent mono- mately must sustain the rapid westerly flow 
tony of surface climatic conditions over wide of higher latitudes. Yet, because of the steadi- 
areas, reflected mainly in small variability of ness of the tropical regime, short-period 
temperature and wind. The opinion prevailed fluctuations of kinetic energy outside the 
that conditions in the free atmosphere were as tropics are not traceable to variations in the 
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Fig. 1. Vertical time section at Eniwetok (11° N, 162° E), September 8—13, 1945 (Greenwich time). A long barb 
represents 10 mph, a short barb 5 mph and a heavy triangular barb so mph. Heavy solid lines are trough- 
lines. 24-hour surface pressure changes are plotted above the surface observations. 


source region of the energy. Poleward trans- 
port of heat energy across the subtropical 
ridgeline is subject only to seasonal and long- 
period variations. 

Gradual accumulation of aerological data, 
however, has revealed that space and time 
variations, especially of the wind flow aloft, 
are surprisingly large. Even the pre-war 
literature carries the annoyed comment that 
the more numerous the observations from the 
tropics, the more complicated tropical ana- 
lysis. The large amount of high-level data 
collected in most parts of the tropics during 
the recent war finally has made it amply clear 
that in many low latitude regions the varia- 
bility of winds aloft equals that of the tem- 
perate zone. 


This is true notably in the upper troposphere. . 


Let us glance at a time cross-section of rawin 
observations taken at Eniwetok in the Marshall 
Island group of the Pacific Ocean (fig. 1), 
a section typical of the sequence of events 
aloft over wide areas of the tropics. In the 


lower atmosphere, up to about 600-mb, the 
steadiness of the flow is undeniable. From 
there on upward fluctuations increase. Extreme 
restlessness rules in the upper troposphere, to 
subside only as we pass across the tropopause 
into the stratosphere. The classical notion 
about the tropics thus is partly right — but it 
gives only half of the story. In the low tropo- 
sphere the course of events is very uniform, 
especially the important intake of heat and 
moisture from the oceans which is dependent 
in part on the surface wind. This, however, 
cannot hold for the poleward discharge of 
heat aloft. At the very least, this heat low must 
vary from longitude to longitude as north and 
south winds alternate along the latitude circles. 

A hemispheric map of the upper troposphere 
(fig. 2) brings this out even more clearly. For 
several months toward the end of the last 
war simultaneous upper-air observations were 
taken throughout the tropical and equatorial 
belts of the northern hemisphere: over Africa, 
the Americas, the Pacific Islands, and south- : 
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Fig. 2. Topography of the 300-mb surface (100’s feet), August 26, 1945, 0600Z. Insert: meridional profile of the 
zonal wind (mps) at 300 and 700 mb computed between longitudes 20° E and 110° E via western hemisphere. 


eastern Asia. These months are eminently 
suited for large-scale analysis. Under auspices 
of the Office of Naval Research, United 
States Navy, the entire period was analyzed 
at the University of Chicago, both at low and 
high levels. The chart for August 26, 1945, 
represents a random selection from this file. 
Even under the favorable circumstances pre- 
vailing in 1945 station coverage was poor in 
some areas, especially in the tropical Atlantic. 
Differential analysis therefore was a standard 
auxiliary analysis tool, together with time 
cross-sections and other techniques employed 
under such conditions. It is not possible to 


vouchsafe for all details of fig. 2. But its 
broadscale features are submitted without 
hesitation as reliable. 

The most outstanding of these features are 
the breakdown of the tropical atmosphere into a 
train of vortices (RIEHL 1948) and the complete 
interlocking of flow between high and low latitudes. 
It would seem more than difficult to draw a 
line separating polar and tropical zones. Every- 
where, troughs of great extent — called 
»extended troughs» by CRESSMAN (1948) — 
reach from pole to equator. It is only to the 
east of these troughlines that heat can flow 
poleward. Thus there is not a general seepage 
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of energy northward, uniform in longitude. 
Heat is injected into the polar zone in few and 
narrow strips of longitude. In part, therefore, 
changes of flow configuration and intensity 
in higher latitudes must depend on the avail- 
ability of low latitude disturbances to form 
extended troughs. We can substantiate this 
point further by calling on an entirely different 
source of evidence. 


Composition of tropical rainfall 


The rate of poleward flow of heat energy 
along extended troughs depends partly on 
the local energy supply in a trough, and the 
energy that can be drawn advectively toward 
it. Since the free atmosphere is a cold source 
with respect to radiation, at least up to 200-mb, 
condensation heat, latent or converted, is the 
only heat energy form available for transport. 
If we should find that the lower atmosphere 
transmits such energy to the upper levels with 
preference along troughlines in low latitudes, 
it would follow that extended troughs— 
whose existence is transitory—are far better 
heat conductors than would be the case if 
advection from large and distant arcas were 
required. 

As the bulk of the low latitude precipitation 
falls in the form of showers, convection is the 
prime mechanism for distributing condensation 
heat vertically. It has been the custom to 
portray convection as occurring predomi- 
nantly with a random spatial distribution, 
limited only by climatic boundaries. This 
pieture follows from the conception of a 
steady tropical circulation. Slight rainfall vari- 
ability should accompany the minute fluctua- 
tions of surface temperature and wind. As is 
well known, this does not hold true, even 
from the seasonal viewpoint. In any given 
arca, excess rainfall may flood the fields in one 
year and drought may strike in the same 
season of the next year. If we determine in 
what manner the individual days contribute 
to the rainfall total, the discrepancy with the 
classical description becomes even more strik- 
ing, even though we hear that in many 
locations showers begin every day at a certain 
hour with clocklike regularity, and even 
though the Hawaiian tongue contains no word 
for “weather”. 

Along with many others who have ventured 
to dispute the classical viewpoint, this writer 


(1945) has maintained that rainfall of the 
tropics is very variable from day to day. At 
any location clear and rainy spells alternate, 
and this alternation occurs in definite relation 
to changes in the upper windfield. Bad 
weather concentrates in narrow zones near 
troughlines. Here we observe high relative 
humidity to great heights, therefore the air 
mass usually labelled as “maritime tropical.” 
Elsewhere, the atmosphere is divided into a 
lower moist and upper dry layer, with a 
boundary between 5,000 and 10,000 feet 
depending on location and season. This 
boundary exists even where the trade wind 
inversion or layer of thermal stability are 
absent. Thus, there are good qualitative indications 
that the upper troposphere receives its latent heat 
supply in large measure in disturbed zones. More- 
over, CRESSMAN (1948) has demonstrated that 
the amplitude of high and low latitude disturb- 
ances and the intensity of the associated bad 
weather increases whenever these disturbances 
join to form extended troughs. 

These statements are susceptible to quantita- 
tive study. The existence of skewness in the 
rainfall distribution at tropical stations on a 
monthly or annual basis is well known. The 
average rainfall is the result of many periods 
with less than average and a few periods with 
heavy precipitation. This skewness becomes 
more pronounced and sometimes extreme if 
we investigate how the daily rainfall contri- 
butes to monthly or annual means. Such studies 
were made recently for rainfall on Puerto Rico 
(RIEHL and SCHACHT 1947) and for Florida?. 
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Fig. 3. Cumulative percentage frequency distribution of 
rainfall at San Juan, Puerto Rico (solid curves) and for 
average of all stations on Puerto Rico (dashed curves): 
a) for dry season, February—March; b) for rainy season, 
August—September. Ordinate: percent of days with 
rain; abscissa: percent of rainfall. Ten years were used to 
compute diagram. Hurricane precipitation is not involved. 

1 H. Remi: “Florida thunderstorms and _ rainfall.” 
J. Meteor. 6, 289—290, letter to editor. 
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Similar results were obtained for the Philippines 
in the course of class instruction. On the island 
of Oahu in the Hawaiian group over 2/3 of 
the annual rainfall stems from about to major 
storms which on the average last 3—4 days 
(RIEHL 1949). As another example, fig. 3 
illustrates conditions at San Juan, Puerto Rico. 
This city enjoys a moderate climate with 
average monthly rainfall ranging from 2.5 to 
7 inches in the course of the year. During both 
dry and rainy season 10 % of all days with rain 
account for 50% of total rainfall, whereas 
so % of all days with rain produce only ro % 
of the precipitation. This holds also for the 
rainfall of Puerto Rico as a whole, including 
many mountain stations. 

It follows that most tropical rainfall occurs 
in consequence of organized, not unorganized, 
convection. Within relatively few and narrow 
zones of disturbed weather does the tropical 
atmosphere actually obtain the largest share of 
that portion of its heat which is derived from 
release of latent heat. These zones, the “‘second- 
ary” disturbances of low latitudes act as 
tropical counterpart of the cyclones of higher 
latitudes in effecting conversion from latent 
to sensible heat. This confirms the hypothesis 
raised at the beginning of this section. The 
bulk of the heat energy transported to high 
latitudes along extended troughs is injected 
into the upper levels within these troughs. 

Our attempt to demonstrate that the tropics 
are a factor in producing changes of the general 
circulation so far has yielded the following: 
the field of motion aloft in low latitudes is 
marked by great restlessness connected with 
the passage of vortex or wave trains; heat 
flows poleward in a few narrow channels, 
variable in longitude and time; the heat 
available for transport is deposited in the upper 
air in a few small areas. These observations 
also imply that it is not possible to arrive at an 
understanding of the general circulation by 
employing only steady state dynamics applied 
to a Statistical mean circulation. This was 
recognized early by Derant (1921) for the 
higher latitudes. RossBy’s studies (cf. Rosszy 
1949) place increasing emphasis on the indi- 
vidual members that compose the mean cir- 
culation. Outside the tropics, the flow of air 
is variable in space and time at all heights. 
It is a peculiarity of the tropics that for the 
most part only the high atmosphere is charac- 


terized by this variability, while the low levels 
are steady. Let us briefly turn to this layer with 
steady circulation. 


The “subcloud” layer 


Perhaps the most perplexing fact about the 
air currents below the bases of the cumuli in 
the tropics is their regular monotone flow. 
Nowhere else is the atmosphere capable of 
maintaining an equilibrium state locally. The 
intake of heat and moisture from the tropical 
oceans, and the heat accumulation in the lowest 
part of the tropical atmosphere actually take 
place in steady broadscale currents, featured 
only by small-scale vertical turbulence. JacoBs 
(1942) has made comprehensive calculations 
of seasonal and geographic distribution of 
energy transfer from ocean to atmosphere. The 
most recent information on the subcloud layer 
itself comes from two Caribbean expeditions 
undertaken by Woods Hole Oceanographic 
Institution. During the first of these voyages 
experiments were made with smoke released 
from ships in order to detect any vertical and 
horizontal eddies that might be present 
(Woopcock and Wyman 1946). Smoke 
plumes were observed to drift downstream 
with lateral distortion of the orientation of the 
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Fig. 4: Vertical structure of the subcloud layer near 

Puerto Rico on April 2, 1946. Airplane ascent was made 

about noon local time. Light dashed line is dry adiabatic 

line, light solid line is line of constant specific humidity. 

Reproduced with permission of Woods Hole Oceano- 
graphic Institution. 
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plumes from the mean wind direction that 
attained over 45°, whereas the calculated 
crosswind varied by 5° at most. In the vertical, 
smoke rose up to 150 m and successive up- 
ward bulges were about 300 m apart. This is 
taken as evidence that the plumes are floating 
in a convection cell pattern that with simpli- 
fication can be patterned on the Bénard cell. 

The second expedition, led by J. Wyman 
and A. H. Woodcock, took ship and aircraft 
measurements of temperature, humidity, tur- 
bulence, and vertical motion from the vicinity 
of the sea surface up to the trade wind in- 
version near Puerto Rico during April 1946. 
A study of the trajectories of the air sampled 
showed it had reached the Antilles on paths 
mainly from the north-central Atlantic (Bun- 
KER, HAURWITZ, Markus, STOMMEL 1949). 


Certain differences compared to the data of 


the “Meteor” expedition therefore should be 
expected. In particular, the large superadiabatic 
lapse rates of the eastern Atlantic should have 
disappeared. This was indeed observed. The 
subcloud layer near Puerto Rico turned out 
to be almost completely mixed (fig. 4) and 
can be regarded as such for all large-scale 
atmospheric problems. In the upper half of the 
layer the potential temperature increases very 
slowly with height (about 0.2°C/roo m). Con- 
versely the mixing ratio decreases. As these 
gradients have inverse effects on the stability, 
the virtual potential temperature also was cal- 
culated. This quantity behaves as the dry 
potential temperature. Thus the upper part 
of the layer has very slight stability with 
respect to dry-adiabatic vertical displacements, 
except where the small-scale horizontal tem- 
perature gradients treated by PRIESTLY and 
SWINBANK (1947) are of importance. In the 
lower portion, vertical gradients could not be 
detected. The temperature at ship’s deck level, 
however, always was lower (about 0.5°C) than 
the surface water temperature. 

The question arises whether the observed 
state of the layer which is essentially steady 
with respect to temperature and moisture is 
due to convection or diffusion. As already 
seen, the Panama smoke experiments favor 
convection as does the temperature difference 
between ocean and deck level. Yet the upward 
increase of virtual potential temperature in the 
upper part of the layer speaks against this 
solution. Moreover, turbulence records taken 
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vapor pressure on April 23, 1946, as obtained on ship 
stationed near Puerto Rico. Abscissa gives time scale 
converted to horizontal distance to indicate eddy sizes. 
The variations of vapor pressure are caused by fluctuations 
of water vapor density as the temperature effect is very 
small. Reproduced with permission of Woods Hole 
Oceanographic Institution. 


in horizontal traverses showed only small 
turbulence gradients, and not even underneath 
the clouds themselves was there an increase 
of turbulence. These latter observations support 
the alternative, that diffusion is the governing 
mechanism. This is also suggested by the fact 
that throughout most of the subcloud layer 
the small-scale fluctuations of temperature and 
humidity mainly are out-of-phase (LANGWELL 
1948) whereas at deck level they are in phase 
(fig. 5). In a convective atmosphere positive 
correlation should prevail generally. 

Fig. 5, however, also indicates presence of 
larger eddies. Superposition of a large eddy 
on more minute fluctuations is easily discern- 
ible. The possibility, if not the likelihood 
remains that vertical motion gradients exist 
even in absence of turbulence gradients. If we 
postulate the presence of patterns of horizontal 
convergence and divergence that have the 
dimension of convection cells, it is still possible 
to seck the “roots” of the cumuli in the sub- 
cloud layer as a whole. 

The Woods Hole observations point to a 
conceivable explanation for the steadiness of 
the windfield in the subcloud layer. Rossby has 
emphasized! that the development of large 
lateral turbulence eddies depends on the ther- 


* Verbal discussions with the writer. The arguments 
based on the relation between vertical and horizontal 
eddies are included with permission of Professor Rossby. 
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mal stratification of the atmosphere. As the 
field of motion always seeks to attain gradient 
wind equilibrium, pressure gradients that can 
sustain the windfield must be able to develop 
aloft. On account of hydrostatic considerations 
this is possible only in the presence of horizon- 
tal temperature gradients that can be intensified 
or decreased. When lapse rates tend toward 
the adiabatic, vertical circulations predominate 
and preclude the formation of horizontal tem- 
perature gradients in the air over the oceans 
that are in excess of the temperature gradient 
of the sea surface. The lateral eddies then must 
recede in importance. This is precisely the 
situation that exists in the trade wind belt. A 
hypothesis can be formulated stating that be- 
cause of lack of thermal stability over nullions 
of square miles, extreme smallness of oceanic 
temperature gradients, and effectiveness of 
vertical turbulence, horizontal circulations will 
be held to a minimum. Flow in the low levels 
should be quasi-uniform, especially under trade 
wind inversions, since they shield the low 
troposphere from any disturbing influence of 
the high levels. 


The cloud layer 


The role of the tropical cumuli in transport- 
ing heat upward has held a place of prominence 
in papers on the general circulation for many 
years. It forms a cornerstone of most theories 
on that subject. Most reasoning is patterned 
to conform to the principles of the “parcel 
method” as prime mechanism of convection. 
The air that has travelled equatorward from 
the polar zones near the surface is thought to 
rise in random convective motion in the 
tropics. These buoyant parcels carry heat even 
to the high troposphere to compensate for 
radiation losses. 

In recent years, dissatisfaction with the parcel 
method has become widespread. In particular, 
the observed temperatures inside cumuli are 
far colder than predicted by this approach. 
Modification of the theory of convection, 
begun with introduction of the “slice method” 
by J. BJERKNES (1938), has been carried for- 
ward recently in a brilliant paper by STOMMEL 
(1947). He explains the observed temperature 
and humidity lapse rates in cumuli with the 
assumption that a buoyant jet entrains air 
from all sides during its ascent. 


Apart from the dynamics of convection, it 
is relevant to ask whether random convective 
heat transfer, as utilized in the older general 
circulation theories, is an effective mechanism 
under any circumstances. Consider a mass of 
air Ms, contained between the top of the sub- 
cloud layer and, say, the 300-mb level. This 
air is to be maintained at an equilibrium mean 
temperature T. If M2, with mean temperature 
T;, is the mass of air added from below in the 
course of one day to offset radiational heat 
losses, the total mass (M) aloft after 24 hours 
must be the sum of M and M2. Mass continuity 
is provided following classical theory by low- 
level inflow from higher latitudes and pole- 
ward outflow aloft. From the law of conserva- 
tion of energy it follows that 


MT SM ATI MET, 


where AT is the heat loss in 24 hours. The 
ratio 


MM = A Te ly, 


where T,— T measures the latent heat of 
condensation realized in the course of one day. 
Let T,—T = 5°C, a large value seldom 
realized even in parcel method computations 
over the pressure interval considered, and let 
AT = 1°C, a moderate amount to judge by 
what is known on radiational cooling in the 
troposphere. Then M2/Mz = 0.2. As the bases 
of tropical cumuli lie near the 950-mb level, 
M = 650g per unit area and M; = 130 g, 
corresponding to a layer of vertical depth of 
130 mb, or roughly twice the total depth of 
the subcloud layer. 

Since vertical mass movements of the magni- 
tude just computed must take place over the 
entire tropics, it is easily anticipated that the 
strength of the average meridional circulation 
will not be capable of maintaining mass con- 
tinuity. In the following section the mean 
meridional circulation will be computed for 
the surface. If the result is valid for the whole 
subcloud layer, the ratio of meridional flow 
observed to that necessary to maintain heat 
balance according to the foregoing is "/ro to 
1/20. It follows that the average equatorward flow 
cannot be a factor of any account in providing a 
mechanism for heat transport from lower to upper 
troposphere. We must fall back on turbulent 
eddy motion. 
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There are two principal choices for the di- 
mension of these eddies. They may have the 
size of convection cells where the vertical 
thermal stratification gives rise to the motion. 
Or they may be “secondary” disturbances 
where the field of vertical motion is controlled 
by the horizontal. The answer has already 
been given. Yet, as the tropical troposphere 
in its entirety is unstable with respect to moist 
adiabatic ascent — apart from trade wind in- 
versions — we might have expected a random 
distribution of cumulonimbi all over the 
tropics, surrounded by clear areas of descend- 
ing air. If this were observed generally, uni- 
form motion might prevail at high levels as it 
does near the ground: Actually the closely 
spaced trade cumuli have their tops between 
6,000 and 10,000 feet. Above this level clouds 
occur rather strictly in ‘organized convection 
zones”, although with “Feinstruktur’’. It is a 


remarkable coincidence that the steadiness of the : 


trade begins to decrease upward rapidly just above 
the tops of the trade cumuli. 

STOMMEL (1947) to date has offered the only 
acceptable explanation for the constantly ob- 
served restraint upon the growth of cumuli in 
spite of existence of positive areas on thermo- 
dynamic charts that often are huge. AUSTIN 
(1948) also has pointed out that because of 
entrainment large density gradients between 
cumuli and their surroundings should not 
occur. An explanation for the calculated en- 
trainment rates has not yet been offered. Yet 
they are evidently of sufficient magnitude to 
prevent growth of random convective cells 
to great heights all over the tropics. The 
average height of the tops of trade cumuli 
outside of areas with strong trade wind in- 
version must be indicative of the rate of con- 
version of heat and kinetic into potential 
energy by entrainment. Because of this con- 
version upward heat transfer in vertical con- 
vection cells does not carry to the upper 
troposphere. This compels the development 
of an entirely different class of disturbances 
—large-scale horizontal eddies in the form of 
vortex trains—which funnel heat aloft in 
zones of organized convection. The vertical 
stability necessary to permit lateral concentra- 
tion of isotherms aloft to support the wind- 
field is produced by the small-scale lateral 
turbulence in the low levels that finds its ex- 
pression in the form of entrainment. 


The average circulation at the surface 


We shall now carry out the 
calculation of the strength of 
the average meridional circula- 
tion at the surface. For this 
purpose, the extensive sum- 


Fig. 6. Percent of circumference of 
globe covered by ocean .between 
latitudes 25° N and 25°S. 6070 
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maries of resultant winds con tained in the 
Atlas of Climatological Charts of the oceans 
are entirely suitable. The region selected 
extends from 25° N to 25°S. About 80% 
of this area is covered by water (fig. 6). 
Since the mean circulation over the continents, 
especially Africa and Australia, resembles that 
over the oceans (cf. BRUNT 1939, pp. 14—15), 
it is reasonable to conclude that the following 
diagrams are valid for the whole equatorial 


belt. 
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Fig. 7. Average meridional wind speed (mps) at the 
surface between latitudes 25° N and 25° S for the area 
given by Fig. 6. 


The strength of the meridional flow (fig. 7) 
which indicates the net non-geostrophic 
mass transport across the latitude circles, is 
I—2 mps or I—2° latitude per day. Air 
crossing the 2sth parallels N and S will require 
2—3 weeks to reach the equator. It is relevant 
to ask whether the curves of fig. 7 represent 
small residual differences between alternating 
regions of northward and southward motion 
around the circumference of the globe on any 
latitude circle. Figs. 8a and 8b show that 
for the most part this is not the case, especially 
in July. In that month, winds with south 
component cover the entire oceanic area al- 
most to 10° N, and the speed of the south 
component exceeds that of the north compo- 
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Fig. 8 a. North and south components of surface wind 
(mps) averaged on each latitude circle for all areas with 
north and south component separately: January. Right 
hand portion of diagram indicates percentage of circum- 
ference on each latitude circle that has south wind. 


Fig. 8 b. Same as Fig. 8 a, for July. 


nent. Flow of mass into the northern hemi- 
sphere amounts to about 1 % of the total mass 
of a hemisphere in 5 days, if the inflow extends 
through a layer soo m deep. 

The boundary between net transport from 
north and south lies near 2° S in January and 
shifts to 15°N in July. This displacement 
furnishes perhaps the fairest measure of the 
seasonal shift of the equatorial axis of the 
general circulation. It amounts to about 20° 
latitude whereas the sun oscillates over twice 
this latitude range. This fact has been taken 
by RossBy (1949) to suggest that the low 
latitude circulation pattern cannot be alto- 
gether thermally determined. Dynamical mod- 
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Fig. 9. Average velocity convergence and divergence 

from 25° N to 25° S for January (solid line) and July 

(dashed line). This diagram was computed by differentia- 
tion of Fig. 7 with respect to latitude. 


ifications must restrict the seasonal circulation 
changes to a smaller latitude range than the 
solar movement. 

Fig. 7 bears out the contention, made 
earlier, that the average meridional flow is an 
insignificant factor for the vertical heat transfer. 
Nevertheless, the computed meridional veloc- 
ity gradients are quite appreciable. This can 
be seen if the curves of fig. 7 are differentiated 
with respect to latitude to obtain the field of 
divergence and vertical motion (fig. 9). The 
magnitude of the divergence is 1076 sec", 
Vertical displacements of 300—500 m per day 
would result at 3 km if the divergence was 
uniform between the ground and this level. 
These values are excessive. It is reasonable to 
suppose that the meridional wind component 
weakens rapidly in the friction layer, a con- 
tention supported by windroses at many sta- 
tions. If the divergence decreases linearly with 
height and becomes zero at 900 mb, the 
vertical motion through the 900-mb surface 
has the order of 100 m/day, a reasonable 
value. 

The principal region of convergence and 
upward motion lies in the northern hemisphere 
in both seasons. It experiences only a slight 
latitudinal shift and a surprisingly small change 
of intensity in the course of the year. A pro- 
nounced secondary convergence zone is evi- 
dent in July but missing in January. It is inter- 
esting to note that surface divergence of con- 
siderable magnitude occurs in both subtropical 
belts during winter. Descent must be taking 
place in the cloud and subcloud layers under- 
neath the trade wind inversion. 

This is confirmed by a glance at the regional 
distribution of surface convergence and diver- 
gence (figs. 10, 11). As should be expected, 
a close correlation exists between the con- 
figuration of these charts and rainfall distribu- 
tion over the oceans. In addition to the broad 
areas of divergence in the subtropics, the 
elongated narrow band of maximum conver- 
gence stands out that marks the equatorial 
convergence zone. Its intensity is almost uni- 
form in both seasons as also noted from fig. 9. 
In the Indian Ocean, a pronounced secondary 
convergence zone extends near latitude 10° S 
in July. There is no trace of a split in January. 
In the Pacific, however, two convergence 
zones are evident in both seasons. Between 
them, the famous dry area of the equatorial 
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Fig. 10. Regional distribution of convergence and divergence (10 ®sec~4) in January. 


Pacific is very well marked. Thus the double The circulation in the upper troposphere 
convergence zone is a local feature of central 

and western Pacific during January. It extends Interest has centered for many years on the 
over half of the equatorial belt in July, while vertical wind distribution over the tropics. 
the other half does not experience it as an Classical theory demands that at a short dis- 
average condition at any time. tance from the equator the easterlies begin to 
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Fig. 11. Same as Figure 10, for July. 
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decrease with height and that a net poleward 
component of motion exists. In view of the 
preceding section, the latter requirement is 
undoubtedly fulfilled. It is not an object of 
this report to enter into the numerous theoret- 
ical arguments that are being advanced to 
explain the observed distribution of zonal 
momentum aloft. Here, we shall inquire 
merely concerning the state of knowledge on 
the actual windfield of the upper troposphere. 
On this point, contradictory statements have 
appeared of late. 

Upper-wind summaries published in the 
course of the years have shown that in some 
locations the lower easterlies give way to upper 
westerlies during all seasons (cf. Fassic 1933, 
STONE 1942). In other regions the easterlies 
attain their greatest strength in the upper 
troposphere at least during part of the year 
(India, Dutch East Indies, Philippines). When 
KUHLBRODT (1928) published his famous dia- 
gram of the vertical wind distribution over 
the Atlantic he stated that westerly winds 


Fig. 12. Vertical variation of the basic zonal flow in the 
Pacific Ocean during September 1945, as determined 
from rawin observations. Diagram from RIEHL (1948). 


replace the easterlies aloft. In the eastern 
Pacific Ocean, wind data at Hawaii show 
predominance of upper westerlies throughout 
the year. 

A recent summary of rawin observations 
over the tropical Pacific during September 
1945 (RIEHL 1948) also indicated upward de- 
crease of the easterlies (fig. 12). HUBERT (1949) 
found strong westerlies at 200 mb in individual 
situations over the equatorial Pacific (Tarawa 
1° N, 174° E; Phoenix Islands 4° S, 172° W). 
On the other hand, computations of a mean 
cross-section along the meridian 80° W led 
Hess (1948) to state: “In low latitudes the 
tropospheric isotherms are practically horizon- 
tal in both winter and summer. This lack of 
sleonoids in the vertical plane would seem 
to render untenable an explanation of the 
maintenance of the trade winds by north-south 
solenoids.” The section of Hess shows easter- 
lies increasing with height in the northern 
hemisphere summer, especially near latitude 
20° N. VUORELA (1948) also has suggested on 
the basis of pilot balloon and radiosonde ob- 
servations taken on a return voyage from 
Europe to South America that there are no 
solenoids to drive the trades. 

This mass of conflicting evidence points out 
clearly that it is very dangerous to draw sweep- 
ing conclusions concerning the tropics as a 
whole from observations taken at a few stations 
and over short periods only. In view of the 
unsteady character of the upper flow, empha- 
sized repeatedly, data such as Vuorela’s hardly 
lend themselves to generalizations. It is neces- 
sary to gain a view on the tropical belt as a 
whole in order to. obtain clarification. On 
account of scarcity of southern hemisphere 
data this is not yet possible. But the series of 
300-mb charts for 1945, mentioned .earlier, 
makes it possible to obtain at least a qualitative 
picture to about 10° N. 

One monthly mean chart—that for August 


1945 


has been chosen for illustration (fig. 
13). In August the zonal circulation in the 
north and the equatorward extent of the polar 
westerlies attain their seasonal minimum. 
Any tendency toward increase of cast wind 
with height in low latitudes should be dis- 
played most prominently. We note on fig. 13 
that even though we are dealing with a 
monthly mean chart, marked contour patterns 
are nevertheless evident. In higher latitudes, a 
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Fig. 13. Topography of the 300-mb surface (100’s feet) for August 1945. This monthly chart was computed from 
31 daily charts. Solid heavy lines are troughlines, dashed heavy lines mark the subtropical ridge. Insert: meridional 
profile of the zonal wind (mps) at 700 and 300 mb for August, 1945, computed for the area given in Fig. 2. 


well-defined westerly current (jet stream) is in 
evidence. Pronounced troughs and ridges 
alternate indicating a wave number of four 
in the polar belt. The outstanding feature ‘of 
the lower latitudes is the array of cellular 
subtropical highs with intermittent troughs, 
in good agreement with the cellular structure 
of the subtropics demanded by V. BJERKNES 
and Cort. (Physikalische Hydrodynamik, 
1933, 679—687). There is no zonal symmetry in 
the sense that all meridional planes show the same 
features. Regions of northerly and southerly tran- 
sport alternate. A general antitrade is not present. 


All major troughs of high latitudes have 
counterparts in the tropics. Even on the 
monthly mean chart the existence of great 
extended troughs is the most impressive 
feature. It is evident, however, that the num- 
ber of troughs in low latitudes exceeds that 
of high latitudes, in part presumably in corres- 
pondence to the much greater circumference 
of the earth. If in addition to the trough posi- 
tions given by fig. 13 we avail ourselves of 
the knowledge that a very persistent upper-air 
trough lies near longitudes 75—80° E during 
summer, we obtain a mean hemispheric wave 
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Fig. 14. Relative topography of the 300 above the 700-mb surface (100’s feet) for August, 1945. ‘“W?’ stands for 
warm and “C” for cold. Insert: meridional profile of relative topography (logarithmic scale) against latitude 
for the area indicated in Fig. 2. 


number of seven or eight for the low latitudes. 
Perhaps the most remarkable asymmetry 
between high and low latitudes occurs in the 
eastern Pacific Ocean. 

It is possible to compute the average zonal 
wind profile from fig. 13, if geostrophic con- 
ditions are assumed. Since the figure repre- 
sents a monthly mean chart, this approxima- 
tion presumably is valid even at low latitudes. 
Many observers have found that except in 
disturbances the deviations from balanced wind 
are smaller in the latitude range 10°—20° than 
previously assumed. Moreover, the geostrophic 
speeds as computed (insert of fig. 13) are in 


good agreement with those obtained from 
rawins in the tropics. Comparison of 300-mb 
and 700-mb shows that the east wind decreases 
with height at all latitudes computed, though it 
is probable from the 300-mb curve that the 
wind direction itself is east at that level at the 
equator. Even if the geostrophic assumption 
is only a rough approximation, the vertical 
change of wind must nevertheless be given 
correctly in the insert. Otherwise it is necessary 
to argue that there are systematic departures 
from geostrophic equilibrium that tend to 
produce supergeostrophic wind at 700-mb 
and subgeostrophic wind at 300-mb with 
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order of magnitude of 5 mps on the monthly 
chart. 

It is possible, however, to dispense entirely 
with wind calculations and merely consider 
the solenoid field. This can be done qualita- 
tively from the relative topography of the 
300-mb surface above the 700-mb_ surface 
(fig. 14). Around the globe there are broad 
areas where the temperature gradient is directed 
northward at all latitudes. Interspersed are 
regions with reversed meridional temperature 
gradient. Fig. 14 perhaps brings out the lack 
of zonal symmetry most clearly. Notably in 
the central Pacific and over Mexico the east 
wind increases with height between latitudes 
20° and 30°. The same holds true over the far 
western Pacific south of latitude 20°. This is 
also confirmed by rawins, especially over the 
Philippines and the Palau Islands group. Thus, 
there is no single meridian that will yield a repre- 
sentative picture of the hemisphere. Fig. 14 is in 
agreement with the cross-section of Hess 
(1948) since along the meridian 80° W the 
warmest air lies near 25° N. But the contour 
configuration over the western Caribbean is 
a local feature only. A little to the east, over 
the Antilles, the east wind must decrease up- 
ward, in accordance with the findings of 
Fassic (1933) and STONE (1942). 

If the contours of fig. 14 are summarized 
for all latitudes, a profile results (insert of fig. 
14) that shows a northward directed temperature 
gradient at all latitudes. This gradient is very 
small in the lowest latitude belt. Perhaps it is 
best to count it as zero there. As just seen, 
however, this does not mean that barotropy 
exists in any particular region. Figs. 13 and 14, 
together with all measured wind summaries 
aloft make it clear that the average flow pattern 
as well as the daily chart is dominated by extended 
troughs and that therefore longitudinal asymmetry 
prevails. Regions of northward and southward 
directed temperature gradient alternate aloft as do 
regions of upper westwind and eastwind. 

Although fig. 14 demonstrates the presence 
of solenoids in low latitudes, it is still necessary 
to determine the energy source that maintains 
the windfield over the tropics. Sometimes it 
has been suggested that this energy may be 
derived from middle latitudes. Calculations of 
the heat energy released through condensation 
and precipitation in organized convection 
zones, however, show that only a negligible 


fraction (5 % or less) of the heat released need 
be converted into kinetic energy in order to 
account for the observed motion. This per- 
centage is very small even in hurricanes 
(about 10%). Thus there is no problem in 
obtaining energy to drive the currents of the 
tropics. The real problem is why there is so 
little motion, why the atmospheric heat engine 
is so inefficient. 


Basic currents of the tropics 


The structure of individual disturbances 
encountered in the variable broad-scale pattern 
must depend on the large-scale features them- 
selves, especially the vertical variation of the 
zonal wind. The westerly current of higher 
latitudes almost always increases with eleva- 
tion. Outside the tropics, therefore, a single 
type of basic current predominates. In con- 
trast, the baroclinity of the easterlies is variable 
as just seen. Moreover, the tropics are subject 
to much larger seasonal changes of basic 
current structure than the higher latitudes. As 
a result, a greater variety of synoptic disturb- 
ances is encountered, and their occurrence 
has a definite seasonal course. In winter 
(generally dry season) the polar westerlies 
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Schematic representation of types of zonal 
currents typical of low latitudes. 
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aloft invade most tropical areas to at least 
latitude 10°. The base of these westerlies slopes 
equatorward with height from the latitude of 
the surface position of the subtropical ridge- 
line. Thus the depth of the trades is variable, 
ranging from a few hundred meters to almost 
the whole depth of the troposphere. Curve ‘c’ 
of fig. 15 illustrates this class of current. 
Curves ‘a’ and ‘b’ show the vertical variation 
of the trade as encountered in the rainy season 
(summer). 

The picture is completed by introduction 
of the low-level “equatorial westerlies”, a 
current that has remained least amenable to 
simple classification and understanding up to 
the present. Originally observed near the 
surface over India, the Far East, and West 
Africa, it was thought to be of southern 
hemisphere origin turning clockwise under 
the law of conservation of momentum as it 
crosses the equator in response to large-scale 
thermal differences between the hemispheres 
(monsoon). While this definition is not in- 
correct for certain areas, it is certainly in- 
complete. Surface westerlies in the equatorial 
zone not only blow toward heated continents, 
but also away from them to their east (RIEHL 
1948). The classical type of equatorial west 
wind should be overlain by easterlies (curve 
‘d’ of fig. 15). This is true in summer over 
India, but it is by no means observed with 
regularity in other regions. For such reasons 
it is best to avoid terminology as “southwest 
monsoon” that has a causal conotation. The 
term “equatorial westerlies”, gives a definition 
that simply refers to the geographic fact that 
these westerlies are separated from the polar 
westerlies by an intervening wedge of easter- 
lies. The interesting fact remains that we 
encounter few areas on earth where west 
winds do not occupy at least a portion of the 
troposphere. 

The geographical arrangement of the four 
currents shown in fig. 15 follows a definite 
pattern. Near the poleward boundary of the 
tropics the trades give way to polar westerlies 
aloft. Low-level westerlies, overlain by easter- 
lies, occur almost exclusively near the equa- 
torial low pressure trough. The heat exchange 
between high and low latitudes of the same 
hemisphere largely depends on the vertical 
structure of the low-latitude windfield, since 
the latter has a definite relation to the meridio- 


nal temperature gradient.! Poleward heat 
transport requires the presence of a northward 
directed temperature gradient, therefore a 
decrease of the easterlies with height. Extended 
troughs form where this temperature field prevails. 
As an example, comparison of fig. 13 and 14 
shows that most of the principal extended 
troughs are situated in areas where the relative 
topography of the 300-mb surface decreases 
northward at all latitudes. In contrast, inter- 
action between high and low latitudes is shut off 
where the temperature gradient reverses. There 
the easterlies increase upward to the tropo- 
pause level. It is typical of these easterlies that 
they are very steady (cf. RIEHL 1948, fig. 9). 
Heat energy realized from condensation pre- 
sumably flows across the equator under these 
circumstances. 

Another characteristic of the basic currents 
deserves mention. This is the striking opposi- 
tion between lower and upper troposphere. 
As indicated in fig. 15, the motion that prevails 
near the ground differs greatly from that of 
the upper troposphere. A finite, generally 
narrow, transition zone extends between these 
regimes. This zone tends to be situated between 
600 and 400 mb, though with many variations 
in individual situations. This observation points 
to the possibility of treating the tropical at- 
mosphere as composed of two distinct layers 
in the vertical, with the soo-mb surface serving 
as dividing line. 

As noted carlier, winds are very steady in 
the low levels while restlessness characterizes 
the high layers. There is an important excep- 
tion to this statement. High-level easterlies 
that increase in intensity with height are very 
steady, as just mentioned, much more so than 
easterlies that weaken with height, or upper 
westerlies. In areas where westerlies occupy 
the low levels and are overlain by casterlies, 
eddy motion generally is strongest near the 
ground. The vertical gradient of wind steadi- 
ness then is reversed. It is noteworthy that a 
reversal of this quantity is always present. 
Pronounced eddies are present either above 
or below soo-mb. They rarely extend through- 
out the entire troposphere, in distinction to 
troughs and ridges of higher latitudes. 


1 This argument is not based on the thermal wind 
equation but on empirical evidence from maps containing 
rawin and temperature observations. The data show that 
the equation is not invalid as far as latitude 10°. 
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Fig. 16. Illustration showing uniform zonal current in 
one layer of troposphere (dashed lines) and equatorial 
shearline and eddy in other layer (solid lines). 


Fig. 16 serves to illustrate this point. The 
pattern shown describes conditions frequently 
observed along the equatorial trough during 
the northern hemisphere summer. Two layers 
are depicted. A broad easterly current occupies 
one of these, while an intense shearline is 
located in the other. Winds are most unsteady 
along the shearline where cyclonic vortices 
develop, and to its south where motion of the 
clockwise eddy usually takes place. In fig. 2, 
conditions in the eastern and central Pacific 
most closely resemble this pattern. The 
astounding fact is that whereas fig. 16 is valid 
almost everywhere along the equatorial trough, 
it must be applied with alternating interpreta- 
tion. In some regions, the shearline is located 
in the upper troposphere and in others at the 
ground (equatorial front or convergence zone). 
In many parts of the equatorial zone the layer 
of steady casterlies adjoins the surface. In 
others it prevails from soo-mb upward. As 
seen from figs. 10 and 11, however, the 
equatorial front zones are very narrow and 
occupy only a small fraction of the total 
tropical arca. 

In view of the foregoing, it is not surprising 
that synoptic studies report gross regional 
differences in equatorial weather patterns. If 
we recognize that a remarkable compensation 
exists between low and high levels, that upper 


and lower troposphere tend to have different 
and often quasi-independent flow patterns, and 
that the relation between these layers varies 
with longitude, it becomes possible to evolve 
a rational account of weather near the equator. 


Conclusion 


This report has attempted to demonstrate 
that the investigator of the general circulation 
cannot proceed like a physicist who turns on 
a well-controlled flame underneath a tank 
with Auid and then studies the motion in the 
tank without further reference to the flame. 
The source of heat for the middle latitude 
westerlies has a very variable character. Trans- 
fer of heat from the oceans surface to the high 
troposphere—the principal layer of trans- 
mittal to the polar zone—is locally and longi- 
tudinally not continuous. This also holds for 
poleward heat emission. Since local and longi- 
tudinal variations of heat exchange are in part 
determined by the tropical circulation itself, 
changes of the flow pattern at the receiving 
end in the temperate zone cannot be under- 
stood entirely without reference to the tropics. 
As yet, no attempt has been made in meteoro- 
logy to evaluate dynamically these poleward 
directed influences. 

Although heat accumulation and emission 
from the tropics is a function of longitude, 
and locally also of time, the question remains 
whether there are fluctuations of these quanti- 
ties in time if we integrate around the entire 
tropical belt. The wartime extension of ob- 
servational material has made possible studies 
of this important problem. Recent research, as 
yet uncompleted, has shown that the question 
can be answered in the affırmative.: The total 
internal energy of middle and upper tropo- 
sphere in low latitudes, therewith the poleward 
heat transport, are subject to marked aperiodic 
fluctuations with a period of roughly two to 
three weeks. This observation makes it clear 
that the tropics must play an important part 
in the production of variations in the intensity 
of the circulation of middle latitudes. 


‘ H. Rrent, T. C. Yen, and N. LASEUR: A study of 
variations of the general circulation (to be published). 
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On the Sensitivity of the Wind Field Relative to 


Pressure Variations 


By SVERRE’PETTERSSEN, U.S. Aır Force 


(Manuscript received 27 October 1949) 


Abstract 


It is assumed that the motion is frictionless and adiabatic but not necessarily horizontal. 
A geostrophic approximation of second order is introduced, and a measure of the sen- 


sitivity of the wind field relative to pressure variations is obtained. 


It is shown that the 


sensitivity depends essentially upon the lateral shear of the geostrophic wind and the 


tangential and orthogonal curvature of the geostrophic streamlines. 


In particular, when 


the geostrophic wind is large, the effect of the orthogonal curvature is predominant. 


1. Notations. — The following symbols will 
e used. 
= (u, v, w) = wind velocity. 
= (4 v, w) = acceleration. 
y, w) = derivative of acceleration. 
= (ug, vz) = geostrophic wind. 
= tangential curvature of geostrophic 
streamlines, positive or negative accord- 
ing as the curvature is cyclonic or anti- 
cyclonic. 
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K, = orthogonal curvature of geostrophic 
streamlines, positive or negative ac- 
cording as the streamlines converge or 
diverge. 

E = mean radius of the figure of the earth. 

g — acceleration of gravity. 

®  — angle of latitude. 

2 = angular speed of the earth’s rotation. 

7 200, f’ = 202 cos @ 

p = atmospheric pressure. 

& = specific volume. 

T = temperature. 

cy = specific heat at constant pressure. 

Ze = height above the reference level of an 
isentropic surface. 

y =the Montgomery acceleration poten- 


tial for an isentropic surface (1937). 
== Ch I ah &Zo 


In particular, if the y-axis is chosen along the 
direction of the horizontal pressure force 
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Furthermore, the following auxiliary nota- 
tions will be introduced 


I dy 
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and may be identified with the components of 
the so-called isallobaric wind, in the meaning 
defined by Brunt and Douctas (1928). 

2. The Equations of Motion. — Choosing the 
x-axis toward the east and the y-axis toward 
the north, the equations of frictionless motion 
may be written 


PERS 
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where X and Y represent the horizontal com- 
ponents of the pressure force. 

Differentiating the above equations with 
respect to time and eliminating # and », we 
obtain 


i + fu — of + wf’ + fiw = X + fY | 


heme ieee AR gay 
Since no further differentiation will be per- 
formed we may leave out of consideration 
certain small terms the magnitudes of which do 
not exceed about one or two percent of the 
major terms. 

Now, f = vf’/E, and f’ = — vf/E. Except 
near the equator and near the poles, fand f’ are 
about equal and of the order of to~4 sec.-1, 
while E is 6.4 + 106 m. In large scale currents, u 
and v are of the order 10 m/sec., w of the order 
Lo mec nd w= 710 © m/sec.2: or less: 
It follows then that the magnitude of the terms 
(fu — fu —f'w) and (fu + ff'w) are about one 


or two percent of those of the terms f?u and 


2y, respectively, and may therefore be omitted. 
J"v, resp y ) 


The differentiated equations of motion, there- 
fore, reduce to 


ii + fu=X+fv| 
| (6) 
Deren | 


We may now choose the system of coordi- 
nates such that the positive y-axis coincides with 
the horizontal pressure force, for the only 
terms in (5) affected by this rotation are those 
involving f’, which have been shown to be 
negligible. With this choice of coordinate axes, 
the horizontal pressure force may be expressed 
in either of the following forms 


ne. fo | 


ox 
(7) 


Now, if the motion is adiabatic, then without 
assuming horizontal motion 


Py ay ay 
ee or" | 
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Substituting from equations (7) and (8) in 
(6), and using the notations defined in section r, 
the differentiated equations of motion may be 
written 


u =Au+Bv+E 


(9) 


where the coefficients have the following 
meaning : 


A=—f(f+V,K); 


D-—/( 


ri, 


v =Cu+Dv-+HR 


BAGGER 


WV, 
von E=f? (V, +1); 
dy 


If the corresponding coefficients are com- 
puted for the isobars in a level surface, additio- 
nal terms appear, expressing the effect of 
the isosteric-isobaric solenoids. Except in 
frontal regions, where the degree of baroclinity 
is appreciable, the difference between. the isen- 
tropic and isobaric coefficients is small, and for 
most purposes the two sets may be used inter- 
changeably. 

From a formal point of view, equations (9) 
are particularly well suited for analyses of the 
oscillations and stability of air currents. In the 
first approximation one may then assume that 
A, B, C, and D (which are second space de- 
rivatives of the pressure force) are constants. 
However, the coefficients E and F (F = v, + 
+ I,, with v, = 0) are essentially first space 
derivatives of the pressure force, and in order 
to obtain solutions to (9) it is necessary to 
introduce assumptions as to the relation betwen 
E, F and u, v. Since any such assumption, 
capable of a clear physical interpretation, is 
difficult to formulate in a manner consistent 
with a complete set of equations, solutions to 
(9) will not be sought. Instead, an assumption 
will be made concerning the acceleration, so 
as to obtain an approximation to the instan- 
taneous relation between the wind in a general 
pressure field. 

3. Geostrophic Approximation of Second Order. 
— It is well known that the assumption of 
geostrophic balance gives a fair approxima- 
tion to the true wind in the large-scale pressure 
systems. It is then assumed that the acceleration 
is small in comparison with the Coriolis 
acceleration, viz., 


|V] <|2QxV| (10) 
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If this is to hold not only at a given moment 
but all the time, it is reasonable to assume that 


IV7|< 4 @OxV)H (11) 


We shall refer to (10) and (11) as geostrophic 
approximations of first and second order, 
respectively. Since neither of the above assump- 
tions can be justified by reference to mathema- 
tical argument, their usefulness must be deter- 
mined by comparison of derived results with 
observation. 

In the following it will be assumed that the 
motion is sufficiently geostrophic for the terms 
containing second derivatives of the velocity 
to be neglected in comparison with the deriva- 
tive of the Coriolis force. 

4. Sensitivity of the Wind Field. — Returning 
now to equation (9) and introducing the second 
order geostrophic approximation, we obtain 
(since B = C) the following expressions for 
the velocity components 


Er DE 
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(12) 
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Now, since the x-axis is chosen along the 
geostrophic wind, v is the cross-isobar com- 
ponent of the actual wind. In balanced motion 
this component vanishes (i.e., CE=AF). How- 
ever, for any value of (CE— AF) different 
from zero, the cross-isobar component and the 
tangential acceleration increases as the denomi- 
nator decreases. Hence, with some justification, 
the quantity 
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(13) 


may be taken as a measure of the sensitivity 
with which the wind reacts relative to pressure 
variations. 

It should be emphasized, however, that (12) 
and (13) represent an approximation, the 


ph 


accuracy of which depends upon the smallness 
of ii and ÿ as compared with uf? and vf?, re- 
spectively. 

It is of interest to note that the above expres- 
sion for the sensitivity of the wind field is 
related (inversely) to the degree of dynamic, 
stability. SOLBERG (1936) considered a zonal 
geostrophic current (i.e. K; = K, — 0) and 
found the degree of dynamic stability expressed 


by (: A il which agrees with (13). VAN 
MIEGHEM (1944) considered a symmetrical 
circular vortex (i. e., K„ = o) and found the 


degree of stability expressed by 


RES ear 
A? 


which (on the same assumptions) differs from 
(13) by the absence of the quadratic term 
ake a 

dy 


= VK; 
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5. Types of Pressure Fields. — Substituting the 
definitions of the coefficients A, B, etc. given 
with equation (9) into (12), we obtain for the 
velocity components 

1 OV, I 
) I 


V, + Ix) { 1 = 

+ (1-458) 4,2 
I I OV. I 

ge —V, K:)(: F ) © Ve K,? 

( Te Éd Sd 


RG 


u = 


(14) 


It will be seen that the deviation of the wind 
from the geostrophic wind ‘depends partly 
upon the isallobaric gradient (here represented 
by I, and I,), and partly upon the lateral shear 
of the geostrophic wind and the tangential 
and orthogonal curvatures of the streamlines 
of the geostrophic wind (i.e. the curvatures 
of the lines p — constant in the barotropic case, 
and y = constant in the baroclinic case). 

(x). Stationary geostrophic wind. — In this 
case I, = I, = 0, and the angle (B) of deflec- 
tion of the actual wind from the geostrophic 
wind is given by 
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PK, 
tan B = — a 
£ 
f— 


dy 


(15) 


As will be seen from Table 1, the magnitude 
of B may be appreciable. At such low geostro- 


Table 1 
Magnitude of the Deviation of the Direction of 
the Actual Wind from that of the Geostrophic 
Wind (in degrees) for Various Values of LE (in 
m/sec.) and K, (in km) in Latitude 45°, com- 
puted from Eq. 15. 


oVg I OVg OVeg 
=-f m = =—f 

oy 2 oy oy 
10 40 - Solo 40 8olro 40 80 

+ Kn} 

400 27 65 76114 45 6517 27 45 
800 14 45 65/7 27 45| 3 14 2 
1 600 45 T4 527) 2 74 


phic velocities as occur near the earth’s surface, 
P is normally less than the frictional deflection 
and is, therefore, difficult to detect. In the 
upper troposphere in middle latitudes where 
geostrophic winds in excess of 50 m/sec. are a 
normal occurrence, very large angles of deflec- 
tion will occur unless the isobars are almost 
strictly parallel. 

It will be seen from (13) that the sensitivity 
decreases as V,2K,2 increases.t Hence, the 
strong winds present in the upper part of the 
troposphere are extremely sensitive relative to 
isobar convergence and divergence. It is often 
observed that appreciable developments in the 
pressure field occur at and below regions of 
such convergence in the upper troposphere. 
Normally, appreciable isobar convergence, or 
divergence, associated with strong winds 
result in pressure variations, with the result 
that the angle B is modified by the isallobaric 
gradient. The figures in Table 1 suffice, 
however, to demonstrate the importance of 
K, as far as the sensitivity of the wind field is 


1 The abnormal case, when either of the two 
parentheses in (13) is negative, will not be considered 
since the wind could not then be predominantly 
geostrophic: 


concerned. Since the acceleration tangential 
to the isobar is equal to fv = fu tan B, it 
follows that regions of appreciable isobar 
convergence or divergence are characterized 
by large accelerations.? j 

(2). Straight and parallel isobars without 
geostrophic shear. — In this case K;= Ky = 
oV,/dy = 0, and (14) reduce to 


PR RE RM (16) 
which is the well-known BRUNT-DoucLAs 
equation (1928), showing that for this type of 
pressure field the actual wind is equal to the 
geostrophic wind supplemented by a vector 
which is proportional to the isallobaric gradient. 
It will be seen that in this simple type of pres- 
sure field S = 1. 

(3). Straight and parallel isobars with geostrophic 
shear. — In this case (14) reduce to 


b= V+ bv =) (— 75) (17) 


It will be seen that the component along the 
isobar is the same as in the foregoing case, 
while the cross-isobar component is larger or 
smaller according as the geostrophic shear is 
anticyclonic (> 0) or cyclonic (< o). It follows 
then that the wind field is more sensitive to pressure 
changes when the geostrophic shear is anticyclonic 
than it is when the shear is cyclonic. 

(4). Curved concentric isobars. — In this case 
K, = 0, and (14) reduce to 


Vik )(s ne ia V;Ki) 3 | 
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and the accelerational components are 
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It will be seen that the motion is more sensitive 
to pressure variations when the geostrophic shear 


and curvature are anticyclonic than when they are 
cyclonic. 


(19) 


Ve K; 


2 SCHERHAG (1948), who places much emphasis on 
the isobar convergence and divergence, interprets 
these features as ccnvergence and divergence of the 
actual motion. 
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The relation between the wind and the 
above-mentioned simple pressure fields is 
shown diagrammatically in fig. 1. In diagram A, 
there is no geostrophic shear and the geostrophic 
departure is determined exclusively by the 
isallobaric gradient, according to Brunt and 
DoucLas equation. In diagram B the geostro- 
phic departure is amplified through anti- 
cyclonic shear, and in C it is reduced through 
cyclonic shear. Diagrams D and E show the 
conditions around a cyclone and an anticyclone 
moving from the left to the right. Here, the 
effects of shear and curvature combine to give 
a large variety of geostrophic departures. It 
will be seen that the winds are rather insensitive 
to pressure changes along the periphery of the 
anticyclone and in the central part of the cy- 
clone, while the reserve is true along the per- 
iphery of the cyclone and near the centre of the 
anticyclone. 

(5). Converging and diverging isobars. — In 
this case the relation between the wind and the 
pressure field is given by (14) and since both u 
and v depend upon all the characteristics of 
the pressure field, the relationship is not 
readily represented in a simple diagram. It 
may be noted, however, that the cross-isobar 
components of the wind are determined pri- 
marily by the isallobaric gradient and the 
orthogonal curvature. These components are, 
however, amplified or reduced according as 
the sensitivity factor is larger or smaller than 
unity. It will be seen from (13) that the effect 
of isobar convergence or divergence (i. e., K;) 
is always to increase the sensitivity factor. 
Particularly in the upper part of the tropo- 
sphere in middle latitudes, where the geo- 
strophic winds are very strong, the sensitivity 
of the wind field may be greatly increased 


[ [ | 


Fig. 1. Illustrating the relation between the actual 
wind (V), the geostrophic wind (J) and the isallo- 
baric wind in relation to isobar shear and curvature. 
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Upper chart: 12-hour tendency of the height 
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even through moderate amount of isobar 
convergence and divergence. 

Figure 2 shows, as an example!, the stream- 
lines of the actual wind in relation to the pres- 
sure force and the isallobaric field. An examina- 
tion of the charts will show that eq. (14) is 
capable of accounting for all essential features 
of the deviation of the actual wind from the 
geostrophic wind. 

The extent to which the sensitivity factor 
(eq. 13) will prove useful as a synoptic indica- 
tion of development in the pressure field 
remains to be determined. It is hoped that a 
synoptic investigation by A. F. GUSTAFSON and 
L.C. Starrett, which will be completed shortly, 
will clarify certain aspects of this problem. 

1 The streamline chart is reproduced with the 
kind permission of Major A. F. Gusrarson, U.S. Air 
Force, from his unpublished report on a detailed analysis 


of the wind and pressure fields over the USA during 
the International Aerological Days April 1948. 
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Duration of Lightning Strokes and Occurrence 
of Multiple Strokes 


By HARALD NORINDER 
Institute of High Tension Research, University of Uppsala 


(Manuscript received 10 December 1949) 


Abstract 


It is shown that the total duration of lightning strokes, the occurrence of multiple 
strokes and the time interval between them are partially influenced when measurements 
are made on direct strokes to tall objects. By recording the radiated electromagnetic 
field of lightning discharges in antenna circuits with cathode ray oscillographs at great 
distances from the lightning paths it was possible to obtain extensive data on the varia- 
tion characteristics of lightning. Special methods were developed in order to analyse 
the lightning discharges with regard to the sequence of multiple strokes and the slow 
moving after-variation. A comparison of the discharge characteristics of lightning in 
Sweden with corresponding results from lower latitudes leads to the assumption that 
the intensity of thunderstorms is greater in low latitudes than in high latitudes. 
This involves researches of a geophysical problem of greatest interest from different 


points of view. 


Introduction 


During the decade before the second great 
war researches were carried out in several coun- 
tries of the electrophysical conditions of 
lightning discharges. In many cases it was in- 
dustrial electrotechnical enterprises which 
were responsible for this research work. Espe- 
cially was this the case in the United States of 
America. As a matter of course the researches 
were to some extent carried out with reference 
to protective applications — hence directed 
on certain practical problems of lightning 
discharges. This did not prevent a more 
general treatment of the problem, which may 
be exemplified by the contents of a compre- 
hensive article by WAGNER and Mc Cann [1], 
entitled Lightning Phenomena. This article, 
which is of great interest from geophysical 
points of view, deals with results up to 1944. 

As a result of the investigations described 
the authors consider certain characteristics of 
lightning that may be said to be quite defi- 


nitely known. On the other hand, others are 
considered to be only partially answered. 

To the first group they consider as solved 
the mechanism of the discharge, its polarity 
and its wave or impulse type. To the same 
group of solved problems they count the total 
duration of lightning discharges, the number 
of discharges in the lightning path and the 
time interval between them. Especially the 
three last mentioned characteristics of light- 
ning discharges are of the greatest importance 
with regard to protective devices against 
lightning strokes. On the other hand, they are 
of the greatest interest from a purely electro- 
physical point of view. 

Against the quoted authors’ opinion that 
the last mentioned characteristics of lightning 
discharges have reached a definite and con- 
clusive solution, I am of the opinion that the 
researches upon which the conclusions are 
founded must not be considered as definite. 
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Fig. 1. Original cathode ray oscillogram of lightning 

discharge taken at the Uppsala Research Laboratory 

twenty years ago. The figure shows predischarges A—B 

and a steep, fronted main discharge with following 
multiple strokes. 


The researches have not been carried out 
during such general conditions that we will 
be able to consider the results as valid for 
general discharges in open country topographi- 
cal conditions. The reasons for this are evident. 
To a great extent the quoted results have 
been obtained through the study of lightning 
strokes to tall structures. In other cases the 
discharge intervals have been determined by 
photographic methods without any considera- 
tion of the influence of the after glowing by 
ionization in the lightning path. In as much as 
the equivalence between the discharge time 
interval as carried out by photographic 
methods, and the corresponding measurement 
by electromagnetic methods has not been 
proved, we must critically consider the former 
method. The weakest point of the results 
hitherto gained must be, as has already been 
pointed out, that the measurements have been 
made on strokes to tall structures. They 


produce very pronounced and local disturb- 
ances of the electromagnetic field from which 
must follow lightning discharge variations 
typical for such a structure but not for open 
field conditions. The above mentioned authors 
themselves remark that the quoted results 
from tall structures are debatable from the 
considered point of view. The tall structure 
exposed in the field of a thunderstorm cloud 
will cause that a greater volume of the cloud 
will be ionized as compared with field con- 
ditions on open and flat ground and hence 
cause a heavier current in the lightning paths. 
This involves both a greater number of 
partial discharges and increased duration. This 
is evident the more as some of the discussed 
measurements have been obtained by using 
discharges to an earth-connected captive bal- 
loon at heights of soo to 800 metres which 
must introduce very pronounced local disturb- 
ances of the electromagnetic field. Similar 
disturbed field conditions even more difficult 
to overlook follow from results obtained by 
measuring lightning strokes to a sharp pointed 
peak of a mountain. Graphs illustrating the 
discussed results will be quoted further on. 

Considering these limitations of the method 
of measuring lightning strokes to tall struc- 
tures it is astonishing that a more safe and 
convenient method to measure lightning dis- 
charges has not been extensively tried. Such a 
method is at hand by using the radiation of the 
electromagnetic field from lightning discharges. 
In this case it is possible to operate with open 
or closed antenna circuit systems — a method 
that has been extensively applied at the Insti- 
tute of High Tension Research. For further 
information about the principles of methods 
see the authors’ article in this journal, vol. 1, 
no. 2, [2] and NoRINDER-DAHILE [3]. 


The General Aspects of the Lightning Dis- 
charges as Investigated with Cathode 
Ray Oscillographs 


Before going into the details of the varia- 
tion features of lightning discharges it will be 
necessary to consider the three characteristic 
stages of the lightning discharges: a) the initial 
— or the predischarges with a total duration 
of up to 20.000 usec.; b) the rapid main and 
partial discharges with a total duration of 
up to 2.000 usec.; c) the slow moving after- 
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variation process with a total duration of 
about 10.000 usec. or occasionally longer. 
Very often we have observed that the total 
duration is shorter than the values mentioned 
above. With the mentioned values at hand 
we have to estimate the total duration of a 
lightning stroke procedure to be of the order 
of up to 0.03 second. It is not easy to record 
with an oscillograph these long time variation 
characteristics of a lightning discharge with 
a satisfactory resolving of the details. During 
our experiments in Sweden, started nearly 
thirty years ago, we used an open antenna 
circuit in combination with a cathode ray 
oscillograph with a rotating drum arrangement 
inside the vacuum tube. At an early stage of 
our researches these records of the variation 
feature of lightning discharges enabled us to 
discover the existence of the predischarges, 
their amplitudes and their duration. From the 
same records it was also possible to verify in 
a more definite manner our earlier discovery 
that the lightning discharges were aperiodic. 
Barlier it had been claimed that they were 
oscillating of a high frequency and sometimes 
as high values as 10° periods were supposed 
to be characteristic of a lightning discharge. 
In fig. 1 we reproduce an example of such an 
old record where the visible predischarges are 
marked between A—B. 


Experimental Considerations when Measuring 
Variations of Lightning Discharges with 
New Arrangements 


We were during these carlier experiments 
constrained to operate the drum in vacuum 
during a certain time period which necessitated 
that several consecutive discharges had to be 
recorded on the same oscillogram. This must 
cause a certain merging of the variation curves 
of the different discharges. This earlier problem 
of analysing the total variation traits of light- 
ning discharges seemed of different reasons to 
be of important value. Hence quite recently a 
special arrangement was introduced in our 
cathode ray oscillographs in order to take up 
these older investigations in a modernized 
way. The cathode ray oscillographs were 
provided with a time base variation sweep of 
an axial length of 20 centimetres. With such 
a linear time circuit, sweeping several times 
above the long time axis, it was possible to 
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Fig. 2. Oscillogram of lightning discharge taken at one 
of the Institute’s field station with recently developed 
recording device for long durations. Predischarges be- 
tween A—B, main discharge at C followed by partial 
discharges from C, to Cy. Long after-variation. 


obtain detailed records covering the whole 
discharge period of lightning. The relay 
arrangement of the author’s cathode ray 
oscillograph construction (ACKERMANN [4]) 
prevented a recording when the cathode ray 
oscillograph was passing the regions of zero. 
This part of the recording is sometimes in- 
terpolated by dotted lines in the following 
curves of lightning discharges. The oscillo- 
graph, when operated with the long time axis, 
was connected to an open antenna circuit and 
the observations were taken with the E method 
which thus resulted in direct records of the 
field force in volts/meter. (See ref. [2]) Some 
results of such investigations are reproduce 

in fig. 2—6. The distances between the light- 
ning paths and the observation stations vary 
between 18 to 40 km, and when possible to 
observe, these distances are given (in km) in 
the figures. A control that local influences did 
not affect the results was introduced by re- 
cording the discharges simultaneously on two 
stations at a distance of 18 km between them. 
This arrangement also made it occasionally 
possible to record the same lightning discharge 
at the two stations. As will be observed later 
on, the amplitude of the electric force of the 
predischarges of a lightning stroke is much 
smaller than the corresponding values of the 
main discharge. Hence it is not always possible 
to obtain a balance value in lightning discharges 
of the amplitudes mentioned such as can be 
recorded in full with a voltage sensibility of the 
oscillographs suitably fitted to one of the two 
amplitudes. Field variation curves drawn from 
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an oscillogram taken with the method men- 
tioned above at one of the field station 
of the Institute, Akerby, at a distance from the 
lightning path of 20 km, is reproduced in fig. 
2. The oscillogram does not show the pre- 
discharges, because they were too faintly 
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Fig. 3. Simultaneous records of lightning discharge from 
two field stations. Typical recorded predischarges before 
main discharge A. No partial discharges are following. 


recorded. The predischarges were thus esti- 
mated to have a duration of 2 msec (milli- 
seconds) marked A—B fig. 2. The main 
and partial discharges and the slow moving 
after-variation are reproduced from the original 
oscillogram. The predischarges are followed 


by the rapid main discharge C and first partial: 
discharge at Cr. On the slow moving after- 


variation superimposed traces of subsequently 


occurring partial discharges are marked as 
C.—C,. The discharges in the lightning path 


which are carrying the heavy currents are 
evidently the ones marked with € and C2 to 
C4. Most frequently it is so that the greatest 
amplitudes are reached by the heavy current 
carrying discharges in the first lightning stroke. 
Sometimes we observe exceptions from this 
rule. See fig. 3, where simultaneous records 
from the two stations at a distance of 18 km 
apart are drawn. The typical predischarges 
with a duration of up to 1.1 msec are visible. 
The lightning is in this case characterized by 
only one main discharge marked with A and 
followed by a slow moving after-variation. 
In fig. 4 oscillograms of three separate 
lightning paths with unipolar variation type 
are reproduced. We observe how variable are 
the durations of the partial discharges or 
multiple strokes. In A: and A; all partial 
discharges are concentrated within a time 
period shorter than a millisecond and in A: 
they extend over more than two milliseconds. 


Oscillograms of two lightning paths ocurring 
within a short time interval and recorded on 
the two stations and characterized by multiple 
strokes are reproduced in fig. 5. How shifting 
the variation type of a lightning may be, is 
shown by a comparison of the simultaneous 
records of the oscillograms of the two stations 
in fig. 6 with the foregoing curves. In this, as 
in the foregoing cases with exception of fig. 3, 
the predischarges were developed with such 
a low amplitude that they failed to be recorded 
within the sensibility of the oscillographs used. 
Hence the time is marked to begin with the 
first main discharge Ax, which is followed by 
6 discharges up to A,, of which the third 
attains the highest amplitude. Later on some 
additional discharges are visible superimposed 
on the slow moving after-variation. 

The typical after-variations of sometimes 
very long duration are caused by rapid disloca- 
tions of charges within those regions of the 
thunderstorm cloud which are taking part in 
the charge transportation by the main and 
the partial discharges. Hence this procedure is 
provoking slow moving after-variations within 
the regions of the thunderstorm cloud which 
are active in the transportation of charges by 
the lightning. If we suppose the transported 
charges to be distributed within homogeneous 
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Fig. 4. Typical unipolar lightning discharges recorded 
from the Institute’s field stations. Limited discharge time. 
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volumes of a sphere, a cylinder or a cone, it is 
rather easy to calculate the field changes and 
their polarity. The above mentioned varia- 
tions of the thunderstorm fields have all been 
measured above open, flat extended surfaces. 
Thus local disturbances in the electric field 
conditions are eliminated. 

If we, on the other hand, consider the corre- 
sponding variations when we are dealing with 
the field conditions in the vicinity of tall 
structures, other field variations are developed. 
As a prominent influence of the intensified 
field conditions in the vicinity of the tall struc- 
tures, we have to expect both a lengthening 
of the time variation of the slow moving after- 
variation in the lightning discharge and also 
an intensification of the field charges. This 
will result in increasing values of the currents 
in main and partial discharges. 


Results from Variations of Magnetic Field as 
Produced by Lightning Discharges 


A typical consequence of such extended 
discharge periods in the lightning paths must 
be that transported electric charges, in the 
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Fig. 5. Two separate lightning discharges recorded at 
short time intervals at two field stations. 
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Fig. 6. Simultaneously recorded lightning discharge at 
two stations. After main discharge A, six partial discharges 
up to A,. Long after-variation at one station, at the 
other station it is hidden by zero band arrangement. 
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Fig. 7. Current variation diagrams of concecutive light- 

ning strokes in same path. Typical diminuation of ampli- 

tudes from first occurring discharge. After ionization 
effects at c and especially d. 


lightning attain very high values of the order 
of 100 coulombs or more. On the other hand, 
the author’s measurements in open field 
conditions by using the electromagnetic frame 
aerial method (ref. [3]) resulted in values of 
up to 10 coulombs and only in some excep- 
tional cases somewhat higher values. Accepting 
the transported high charge values, e.g. of the 
order of 100 coulombs, as generally valid also 
in lightning discharges under open field condi- 
tions one is led to field force values within the 
charge delivering thunderstorm cloud of 
quite another order of magnitude than the 
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one generally accepted. This discrepancy is 
not eliminated by the absurd and against all 
experience contradictory explanation that the 
surface of the thunderstorm cloud which 
should be necessary to maintain a lightning 
discharge, must be 100 square kilometres. 

The above mentioned tendency of the 
maximal current values of successive strokes 
in the same channel to decrease in most cases 
can easily be explained. The successive trans- 
portations of charges in the multiple strokes 
must diminish the available charge within the 
regions of the cloud which are active in the 
discharges. We have another way to check 
this typical decreasing tendency of the current 
intensity in lightning discharges. An extensive 
investigation of the variation of the mag- 
netic field as caused by vertical lightning dis- 
charges has been carried out in Sweden (ref. 
[3])-. Under certain conditions these researches 
have allowed a calculation of the current 
variation of successive discharges in the same 
lightning channel. The decreasing amplitude 
values are illustrated by fig. 7. In the vertical 
lightning from which the 4 discharges (a—d) 
have been oscillographed, they are evidently 
passing consecutively in the same path. Of 
special interest is that consecutive discharges 
in the lightning channel very often possess 
similar variation forms. From an electrophysi- 
cal point of view this is quite easy to explain. 
The reason is that the transportation of the 
charges takes place in the same channel and 
hence is fed from the same ionized and active 
region of charges of the thunderstorm clouds. 
Worthy of note is the retardation of the steep 
front especially in d, where a weaker current 
amplitude precedes the rapidly increasing 
current front. Evidently this must be explained 
as being caused by a retention of the ionization 
from the foregoing discharge in the same 
channel. Very characteristic is the gradual 
diminuition of the amplitude of the consecu- 
tive current impulses and hence a diminuition 
of the transported charge. This is easily 
observed in the curves a — e.of fig. 8, which 
represent the consecutive current discharges 
in another vertical lightning path. 

Of special interest in this connection is the 
distribution of the total durations of main and 
partial discharges which constitute vertical 
lightning paths. Such an analysis has been 
carried out for 106 individual discharges in 


the above mentioned earlier investigation at 
the Institute of High Tension Research (ref. [3] 
and fig. 9). The real duration values are in this 
case insignificantly longer, depending on the fact 
that weaker parts of the recorded field force 
in the oscillographs are to some extent hidden 
by the zero-band of the instruments (ref. [4]). 
The most frequently occuring duration values 
are from 100 to 200 psec. This is in full accord- 
ance with duration values that have been 
measured by quite another method (see fig. 
2—6). In this case atmospherics were taken 
up by open antenna circuits. The duration 
of the rapid main and partial discharges was 
of the same order of magnitude in fig. 9. 
Duration values of up to 1000 psec. are, as 
may be seen from fig. 9, very rare. It is of 
importance to keep the quoted duration values 
in mind when we have to deal with the prob- 
lem of the total duration of lightning strokes. 
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Fig. 8. Current variation diagrams of concecutive light- 
ning strokes in same path. Typical diminuation of cur- 
rent amplitudes from first occurring discharge. 
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Investigations of Radiated Electromagnetic 
Field of Lightning as Necessary for Study of 
Variability of Discharges 


The very pronounced variableness, which 
appears to be characteristic of lightning 
discharges, requires special precautions if we 
are trying to carry out a reliable statistical 
investigation. This must be considered with 
regard as well to the duration of the discharges, 
as to the intervals between them and the total 
duration of lightning strokes. Evidently the 
primary condition must be to analyse by 
cathode ray oscillographic recording methods 
a sufficiently great number of individual 
lightning strokes. As we shall observe further 
on it is rather an easy task to obtain extended 
observations, especially if the researches are 
extended to include lightning discharges 
between clouds, and thus not only limited 
to vertical ones against the earth. In this con- 
nection we are able to use extensive experiences 
from another recent investigation not yet 
published. This one deals with thousands of 
oscillographic records of atmospherics taken 
during the thunderstorm season of 1949. The 
atmospherics were recorded simultaneously in 
Sweden with the aid of two stations at a 
distance of 576 kilometres apart. From this 
extended accumulation of recorded lightning 
discharges we are able to conclude that there 
is no pronounced difference in character, 
either with regard to duration of individua 
consecutive discharges in the lightning strokes, 
or with respect to time intervals between the 
discharges. 

Some examples of typical consecutive light- 
ning discharges recorded from separate light- 
ning paths at one of the two mentioned sta- 
tions at a distance of 576 km apart, are re- 
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Fig. 9. Duration of currents of individual discharges of 106 
vertical lightning discharges recorded at the Institute. 
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Fig. 10. Consecutive discharges in distant lightning 


channels recorded at the Institute’s field station. 


produced in fig. 10—11. In these cases the 
lightning discharges were located at a con- 
siderable distance from the recording station. 
This must to some extent have influenced the 
records by reflection effects from ionospheric 
layers. Such special influences will no tbe 
discussed in this article. From fig. Io—1II, 
follows that consecutive discharges in the 
paths of the examples chosen do not show 
any appreciable differences in their variation 
character. 

We therefore decided to carry out, during 
the latter part of the recent thunderstorm 
scason of 1949 in Sweden, an extended analysis 
of the number of strokes, their total duration 
and the time interval between multiple dis- 
charges in the lightning strokes. For this pur- 
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Fig. 11. Consecutive discharges in distant lightning 
channels recorded at the Institute’s field station. 


pose we chose a certain day, Oct. 20th 1949, 
when we had reason to expect only two re- 
gions of lightning storms in Sweden. They 
were situated at such great distances apart that 
the measurements did not become aftected 
by local disturbances in the vicinity of the 
thunderstorm centres. The map of fig. 12 
illustrates these centres where the filled circles 
represent damages by strokes to earth-connect- 
ed structures, such as buildings, telephone 
lines or electric lines. 

In order to measure the lightning discharges 
from the centres mentioned, an open antenna 
circuit in combination with a cathode ray 
oscillograph was used [2]. The method allowed 
records of the rapid variations of the vertical 
field force as produced by the lightning strokes 
from the thunderstorm centres. If the time- 
base circuit of the oscillograph was adjusted 
to a suitable time scale we could thus record 
that each rapid discharge in a lightning stroke 


—thus the consecutive multiple strokes — 
appeared as rapid vertical deviations from the 
zero line of the oscillograph. These vertical 
deviations were followed by traces of a damped 
after-variation with a gradually diminishing 
amplitude. The direction of decrease of this 
after-variation enabled us to determine the 
number of multiple strokes and the sequence 
of them. The records are of a similar type as 
the ones reproduced in fig. to—11. From the 
above mentioned investigation of records of 
lightning discharges, simultaneously made at 
two stations far apart, followed that all light- 
ning discharges of appreciable intensity passing 
in the channel of a lightning stroke, can be 
recorded provided the amplification of the 
oscillograph is increased to a sufficiently high 
level. 

Hence, the described method should result 
in an extended analysis of the number of 
discharges in lightning strokes. It was also 
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Fig. 12. Map of distant lightning discharges. Filled 
circles indicate damages of strokes to ground. 
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possible to analyse the time intervals between 
the multiple strokes and the total duration of 
lightning discharges. The total amount of 
separate lightning paths, recorded in the 
investigation, was 257. Of these only $7 
consisted of one discharge in the path, 200 
carried two or several discharges of which 
it was possible to measure 686. With the aid 
of these it was possible to determine 629 
intervals. 

The collected and analysed records allowed 
a calculation of the percentage number of 
multiple strokes reproduced in fig. 13. Evi- 
dently, the number carrying only one light- 
ning discharge in the channel must in such a 
distribution become somewhat preponderant. 
In this connection it will be of special interest 
to compare the distribution with the number 
of strokes in 138 vertical lightning paths see 
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Percentage distribution of the number 
of components in lightning strokes 


Fig. 13. Percentage distribution of number of compo- 
nents recorded from 257 separate lightning paths. 
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Fig. 14. Percentage distribution of number of compo- 
nents of 138 vertical lightning paths. 
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Fig. 15. Percentage distribution of total duration of 250 
lightning strokes. 
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Percentage distribution of time intervals between beginnings 


of successive components of multiple lightning strokes. 


Fig. 16. Distribution of time intervals of successive 
components of multiple lightning strokes. 


fig. 14. This distribution has been taken from 
an earlier investigation (ref. [3]). When we 
are considering the shifting variation features 
of lightning discharges in total we must 
accept that the distributions in the two cases 
show a surprisingly good conformity. 
Multiple strokes in excess of four are less 
frequent than the range on four. This is easily 
accounted for, the more so as the charges of a 
thunderstorm cloud must be distinguished by a 
certain limitation. The active charges within 
the cloud transported by a lightning, can in 
some cases be diminished to an insignificent 
value by merely one or a few discharges. 
From this point of view it will indeed be of 
a special interest to compare the number of 
discharges in different types of thunderstorms. 
Especially a comparison of the conditions 
within thunderstorms in the temperate zones 
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Percentage distribution curves of the number of components in lightning strokes 


1 Wagner.Mc Cann & Beck 9 strokes to tall objects 

2 = S = 37 photographs of strokes to open ground 
3 i 5 40 arrester surge currents 

4 Me Eachron 184 protector tube surges 

5 Schonland 65 photographs of sfrokes to open ground 
5 Stekolnikov & Valeev 25 induced surges 

7 Norinder 257 strokes 


Fig. 17. Distribution of number of components in light- 
ning strokes. 


with those in the subtropical zones should be 
of great interest. 

The percentage distribution of the total du- 
ration of the multiple discharges, carried out 
in a logarithmic scale, are reproduced in fig. 
15. The most frequently occuring values are 
varying between 20—40 msec. A percentage 
comparison, of the durations of the intervals 
between multiple strokes in the same channel 
is reproduced in fig. 16. Evidently the dura- 
tion of the intervals is evenly distributed on a 
logarithmic cale. 


Comparison with Results concerning Duration 
and Variability of Lightning Discharges in 
Low Latitude Countries 


It is of special interest to compare the re- 
sults obtained by the investigation of the above 
mentioned 257 lightning strokes with values 
obtained by other methods and under different 
climatic conditions. In this respect it is simple 
to try a comparison with the data compiled 
by Wacner and Mc Cann [1]. There are 
several reasons for such a comparison. The 
compiled values have been gathered from ob- 
servations from different countries in the 
world, all of them situated at much lower 
latitudes than Sweden. 

In order to carry out such a comparison the 
curves in the paper by WAGNER and Mc Cann 
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dealing with the number of multiple strokes, 
the percentage distribution of time intervals 
between strokes and the total durations have 
been redrawn and the recent swedish results 
have been entered in the redrawn curves. 
In the original curve systems the authors 
included some measurements from Sweden 
taken from my first observations of magnetic 
field variations from lightning strokes. In this 
case a differential method was used, which 
did not allow of a record of the slow moving 
after-variations of the current pulses. This 
outcome is quite contrary to the integral 
methods introduced in a later publication ref. 
[3]. The quoted older data do not quite 
correctly reproduce the time intervals be- 
tween subsequent current pulses and hence 
they have been excluded from the new curve 
system of time intervals. 

In fig. 17 the percentage number of com- 
ponents in lightning strokes is reproduced. The 
number of strokes from the Swedish material 
does not result in a distribution curve with 
such high values as are valid for strokes to tall 
objects. The curve obtained from our great 
number of observations represents a typical 
average picture and an extension of the obser- 
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Percentage distribution curves of time intervals between the 
beginnings of successive components of multiple lightning 
strokes. 
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Fig. 18. Time intervals between successive components 
of multiple strokes. 
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vations will probably not alter this general 
character. 

In fig. 18 is reproduced the percentage dis- 
tribution curves of time intervals between the 
beginning of successive components. In this 
distribution the Swedish values [4] are of a 
pronounced shorter time variation type, and 
mostly approach to the curves 1 and 2. The 
values of 3, taken in South Africa, have 
pronounced time intervals of longer durations. 

In fig. 19 the percentage distribution of 
durations of lightning strokes is reproduced. 
Also in this case there is a pronounced differ- 
ence between the durations as measured in 
Sweden compared with those in countries 
situated at lower latitudes. 

There is no doubt that the influence of the 
tall objects upon the results cannot be neglected 
when the lightning mechanism and especially 
the time variation character of lightning 
discharges are considered. The differences in 
the curves of fig. 19 are pronounced also when 
we are dealing with measurements where 
strokes to open ground are analysed. These 
differences must have another explanation. I 
am inclined to ascribe the differences to a 
certain extent to latitude effects, as indicated 
previously. I have not personally been in 
a position to investigate lightning discharges 
in countries situated at a more southerly lati- 
tude, but I have had many opportunities to 
observe thunderstorms visually in southern 
countries, both in Europe, Africa and North 
America. My personal impression is that the 
intensity of such thunderstorms at lower 
latitudes in most cases is at a much higher 
level than in Sweden. 


Characteristics of Lightning Discharges as 
Influenced by Latitude Differences 
of Thunderstorm Intensity 


There are several groups of observations 
which in an evident way manifest that the 
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Fig. 19. Distribution of total duration of lightning 
strokes. 


thermodynamic intensity within thunder- 
storm clouds very often must be remarkably 
higher in low latitudes. From this it must 
follow that the amount of energy transformed 
into electic charges within the thunderstorms 
also must attain a higher intensity level. There 
is no doubt that this quality to some extent 
must guide to the physical explanation as to 
why both the lightning current intensity, the 
length of the intervals between multiple 
strokes and the total duration of the lightning 
strokes will sometimes result in marked 
differences when thunderstorms in higher and 
lower latitudes are compared. Obviously it 
should be of special interest to undertake a 
thorough investigation of this problem. The 
problem is of importance with regard to the 
propagation and transformation of atmo- 
spherics at long distances. Another aspect 
must be to what extent lightning currents to 
the earth are contributing to the maintenance 
of the negative charge of the earth. 
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Atmospheric Current 


By P. QUENEY, University of Paris 


(Manuscript received 25 November 1949) 


Abstract 


By a proper choice of coordinates and variables, a simple and quasi-symmetrical system of 
differential equations is obtained for the infinitesimal adiabatic perturbations of the most general 
continuous zonal air-flow. The coordinate surfaces are those of the orthogonal system built up 
with the isentropic surfaces and the meridians, and the variables are the displacements of the 
fluid in the directions of the coordinate lines, and the local disturbance of the pressure. From 
these perturbation equations a differential equation for the pressure disturbance alone is derived 
in the case of a wave, and it is then shown that this wave equation can be always greatly simpli- 
fied if account is taken of the order of magnitude of the various parameters. Simplified expres- 
sions are also given for the components of the vibration (fluid displacement). More particular 
equations and formulae are then derived for short-period, middle-period and long-period waves 
successively, this classification being made according to the value of the orbital frequency com- 
pared to the coefficient of hydrostatic stability and Coriolis parameter, and special cases are 
also considered where the wave equation reduces to one with constant coefficients (elastic, 
gravitational, Rossby’s waves). Finally some incorrect recently published results are discussed. 


I. Definition and meteorological significance 
of the perturbation problem 


It is well known that when a fluid motion 
may be considered as a disturbed motion, that is, 
resulting of an infinitesimal perturbation 
superimposed upon a given undisturbed motion, 
its differential equations can be linearized, in 
other words replaced with a system of linear 
equations, the so-called perturbation equations, 
the resolution of which is of course easier than 
that of the motion equations themselves, at 
least when the undisturbed motion is a simple 
one. 


However, when considered from this view- 
point, the perturbation problem (resolution 
of the perturbation equations under given 
conditions) might look to be one of only 
mathematical interest, and the large place 
which it takes in modern meteorology would 
not be justified. It seems therefore necessary 


to give here a brief discussion about this 
topic: what is the practical significance of the 
perturbation problem in meteorology? 
During the past two centuries, the perturba- 
tion method was used by several of the greatest 
mathematicians or physicists when they built 
the now classical theories of such phenomena 
as sound-waves (NEWTON, LAPLACE, POISSON, 
STORES, RAYLEIGH), occan’s tides (LAPLACE), 
gravitational water-waves (LAGRANGE, CAUCHY, 
POISSON, STOKES, GREEN, AIRY, KELVIN), ca- 
pillary-gravitational waves (Kervin). Nobody 
can deny the great success of these theories 
for the explanation of important observed 
facts, but it must also be noted that most of 
them only concern very simple problems: 
the fluid in the undisturbed state is stagnant 
or in uniform motion, and possesses no inter- 
nal stratification, so that the coefficients of the 
perturbation equations are generally constant, 
and the perturbations are essentially stable. 
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The first example of unstable waves was 
given by HELMHOLTZ 1868 when he considered 
the oscillations of the system formed by two 
homogeneous currents moving side by side 
with uniform but different velocities (shearing- 
gravitational waves). In this case, the shear at 
the common boundary is the only factor of 
instability, and the only stabilizing factor is 
gravity, which operates mairly on long waves. 
From this it follows that all the waves are 
unstable the wave-length of which is smaller 
than a critical value, whence HELMHOLTZ 
concluded that a turbulent layer of transition 
must develop between the two currents, the 
width of this layer being comparable to the 
critical wave-length (in modern terminology, 
the system is said to be “dynamically unstable”, 
or to possess a “dynamic instability”). HELM- 
HOLTZ applied his theory to stationary billow- 
clouds, by assuming that the surface of dis- 
continuity is the upper boundary of an at- 
mospheric inversion. 

In the case of air over water, this problem 
was generalized by Kervin 1871, by introduc- 
ing capillarity .as an additional stabilizing 
factor, and this author made an application 
to the theory of the formation of water-waves 
by wind. Since capillarity mainly operates on 
short waves, and gravity on long waves, it is 
found that the system is dynamically unstable 
only if the wind exceeds some critical value, 
and if the wind increases from this value, 
there is a limited spectral-interval of unstable 
waves, the width of which increases from zero. 
We have here a simple example of an “in- 
stability criterium’? (condition for the un- 
disturbed motion to be dynamically unstable), 
which may be physically interpreted as the 
condition for the motion to change from 
laminar to turbulent type in the vicinity of 
the surface of water. 

If the theory be correct, the critical velocity 
would be that of the slightest wind which is 
able to raise waves on a smooth water-surface, 
and the theory would also give the wave- 
length of these first-appearing waves. As a 
matter of fact, it is actually observed that 
there is a definite wind-velocity at which the 
surface of water begins to be ruffled, but 
this velocity is much smaller than Krıvin’s 
critical value (about six times smaller). The 
explanation of this disagreement, as pointed 
out by JEFFREYS 1925, is probably that a turbu- 


lent air-layer begins to form on the lee-side 
of cach wave at a wind-velocity which is much 
smaller than the critical value, thus increasing 
the effect ofwind-shear as a factor of instability. 
But there is still no complete theory of wind- 
waves. 

It has been noticed for a long time that 
many of the atmospheric perturbations have 
the same character as gravitational or shearing- 
gravitational waves, but this analogy remained 
a purely formal one as long as the general 
structure of the atmosphere was not known. 

The earliest theory concerning atmospheric 
waves was perhaps HELMHOLTZ’S already 
mentioned solution. of the billow-clouds prob- 
lem, and a quite similar theory is that pro- 
posed by KÜTINER 1939 for mountain lee- 
waves (in both cases an inversion boundary 
is assumed to exist). 

But by far the most interesting surface-wave 
theory applicable to the atmosphere appeared 
after J. BJERKNES and his co-workers dis- 
covered the atmospheric fronts in 1919, and 
were able to describe the cyclones of Northern. 
Europe as wavy disturbances growing and 
progressing along a quasi-permanent tempera- 
ture-discontinuity, the so-called polar-front 
surface. This interpretation could not fail to 
suggest a close analogy between cyclone- 
waves and wind-waves for instance, and 
therefore it was natural to consider the theo- 
retical problem of the oscillations of the system 
formed by two air currents having uniform 
but different temperatures and velocities, 
separated by the polar-front surface and 
bounded downwards by a horizontal surface 
at sea level. Unfortunately this problem is 
much more difficult than Kervın’s wind- 
wave problem, specially because of the slope 
of the polar-front surface, which cannot be 
neglected at the scale of cyclone-waves. 
Tentative approximate solutions were given 
by SOLBERG 1928 and then by V. BJERKNES 
and SOLBERG 1933, and the general results 
were quite similar to KELVIN’S ones con- 
cerning wind-waves: owing to the fact that 
there are two stabilizing factors, namely 
the gravity, operating now mainly on the 
shortest waves (the very short waves, of the 
scale of billow-clouds, are disregarded in this 
problem), and the earth’s rotation, which on 
the contrary mainly operates on the longest 
waves, the system is dynamically stable as 


ADIABATIC PERTURBATION EQUATIONS FOR A ZONAL CURRENT 37 


long as the wind shear is smaller than a critical 
value, and when this value is just exceeded a 
dominant. unstable wave appears, which may 
be interpreted as the cyclone wave. 

This remarkable result, following the great 
push of enthusiasm created by J. BJERKNES’ 
discovery, largely contributed to the success 
of the Norwegian meteorological school. 
However, when considered objectively, this 
result is in reality more a disapproval than a 
justification of J. BJERKNES’ fundamental idea 
that cyclones are the result of the conflict be- 
tween two air masses of different temperatures: 
in V. BJERKNES-SOLBERG’S waves, a temperature 
difference is indeed a stabilizing factor, the 
only factor of instability being the wind shear. 
On the other hand, it very soon appeared 
that the polar-front surface is not a continuous 
surface, that it never exists out of cyclones, 
and that many extra-tropical cyclones, for 
instance those of the mediterranean area, do 
not show any front at their early stage. It was 
therefore obvious, in the years preceding the 
last world war, that the problem of cyclone 
formation had to be completely reexamined. 

In presence of this failure, the modern 
meteorologists were forced to abandon the 
idea that atmospheric perturbations are ne- 
cessarily related to discontinuities, and they 
were led to consider, as factors of stability 
and instability, those arising from the internal 
stratification and wind-shear of the atmospheric 
layers. The simplest cases being obviously 
those where there is no internal wind-shear, 
one must expect to find here the first successful 
attempts of this recent stage of theoretical 
meteorology. Let us mention specially the 
beautiful results arrived at by LAMB 1932, 
TAYLOR 1936 and PEKERIS 1937, concerning 
the atmospheric tides, RossBy’s well-known 
theory of long waves in the westerlies, and 
also, as provisional theories, HAURWITZ’S results 
about billow-clouds, and the recent researches 
on mountain waves. 

But it is clear that, if the assumption of no 
shear is maintained, there can be no hope to 
arrive at an even approximate theory for 
atmospheric perturbations of the scale of 
cyclones, since these perturbations are always 
associated with important horizontal tempera- 
ture-gradients, and these in their turn are 
associated with vertical wind-shears. This 
consideration explains why the cyclone prob- 


lem has become practically identified with 
the problem of unstable waves in a baroclinic 
atmospheric current, and why so much 
importance is attached in present days to the 
derivation of instability criteria for such 
currents. Unfortunately, the results obtained 
up to now in this connection are rather dis- 
cordant: according to a first group of meteoro- 
logists (SOLBERG 1936, HOEILAND 1939 and 
1941, FJOERTOFT 1941 and 1946, VAN MIEGHEM 
1944 and 1946), it is possible to derive, at 
every point of a baroclinic air flow, an in- 
stability criterium which only involves the 
local values of hydrodynamical parameters, 
while according to others (JAW 1946, CHARNEY 
1947, EADY 1949), on the contrary, the bound- 
ary conditions play a dominant röle, and for 
the current values of the wind-shear all the 
short waves are unstable as far as the vertical 
acceleration is neglected (the very short 
waves are stable, but they are unimportant 
in the problem), in other words the system 
is practically always dynamically-unstable. 
(In order to reconcile this result with the 
existence of a dominant wave-length in the 
turbulent motion, Eady imagines to introduce 
the postulate that, among the numerous un- 
stable waves, it is the wave of maximum 
growth-rate which must predominate.) 
What may be now the conclusion of this 
brief historical survey? The failure of the 
perturbation method when applied to shearing 
or baroclinic systems is probably to be attri- 
buted, not to the method itself, but well to 
the incorrectness of the assumptions made by 
the authors of theories, concerning the hydro- 
dynamical conditions imposed to the motion: 
in the wind-wave problem, for instance, 
the motion must be taken as turbulent in- 
stead of laminar as in TayLor’s theory, and 
in all of the recent researches about baroclinic 
currents, the simplifying assumptions made in 
course of the computations are generally not 
permissible, as we shall prove it for some of 
them at the end of this paper. It is therefore 
to be expected that, once the perturbation 
problem is expressed in correct analytical 
form, and the mathematical derivation are 
correctly made, the same success will be 
obtained for baroclinic as for stagnant fluids. 
However it is obvious that the perturbation 
method will be always quite unable to explain 
and forecast the growth of perturbations inside 
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the turbulent regions which result from dyna- 
mic instability. This is in fact a very serious 
limitation of the meteorological significance 
of the perturbation problem, and it may be 
that most of cyclonic perturbations, for in- 
stance, have nothing to do with dynamic in- 
stability (this is strongly suggested by J. 
BJErKNEs’ original theory). However there is 
certainly a great deal of meteorological features 
such as long waves, mountain waves, main- 
tainance of turbulence (cf. RossBy 1947), etc., 
which will be largely explained by using 


perturbation method. 


2. General scope of the present paper 


In the preceding section we have mentioned 
that, because of the great difficulties en- 
countered in the mathematical derivations, 
most of the authors of recent researches about 
the perturbations of baroclinic currents have 
been led to introduce several simplifying 
assumptions in the course of the computations, 
and that these assumptions are generally not 
permissible. Since the most interesting prob- 
lems are probably those related to zonal 
currents, we have found it important to give a 
system of perturbation equations for the most 
general continuous zonal stream, and to show 
that these equations can be practically always 
written in a relatively simple form if proper 
variables and coordinates are used, when the 
parameters are taken with the usual order of 
magnitude observed in the lower (or meteoro- 
logical) atmosphere (up to 30 km, say). This 
is the main scope of this paper. In the last sec- 
tions, we also give some applications to simple 
cases, and a brief discussion about some in- 
correct recently published results. 


3. Definition and general properties of 
symbols 


a) Aj, abridged notation for 24/0; 
D/Dt, individual derivative with respect to 
time f 
A, value of A in undisturbed motion 
OA, local disturbance of A, that is, variation 
of A, from undisturbed to disturbed 
motion, at fixed time and point 
DA, individual disturbance of A, that is, 


variation of A at fixed time and fluid- 
element 


the computation rules for differentials 
hold for the operators d and D separately) 
(A), order of magnitude of A 
AB (or B— A); A>B (or B < A); 
A = B (or B= A), means that (A) = (B); 
(A) > (B); (A) > (B), respectively 
A + B (or B + A), means that A—B<A 
b) M, position vector of the current point M 
defined by the space coordinates 44, 9» 43 
(DM is the displacement of the fluid- 
element from undisturbed to disturbed 
motion) 
Dx,, Dx, Dx, components of DM 
along x, X, X3, directions of the three 
coordinate lines 
9/9x,, d/0x,, 9/9x,, partial derivatives in 
Che dite HO avy. Wa, Ga 
DA = A,, Dgı + Agg Dq2 + Ag3 Da; = 
=) 4, Dé + A, De, PASS 
DA DA CHA (32) 
2 (94/9) = (0/9) (04); D(DA/Di) ve 
3-3 


= (D/Dt) (DA) 

c) i, imaginary unit; log A, Neperian logarithm 
of A; A-B, salar product; A x B, vectorial 
product 

d) L, r, h, cylindrical absolute coordinates 
relative to earth’s axis: L, absolute longitude; 
r, distance to earth’s axis; h, distance above 
equatorial plane 


x, y, 2, axes pointing eastwards, northwards, 
upwards, respectively 
x, Y, Z, direct system of rectangular axes 
such that Y is directed northwards parallel 
to isentropic surface 
x, unit vector of x-direction 
e) g, relative gravity (that is, ordinary gravity) ; 
E, potential of absolute gravity; p, pressure; 
q, density; T, absolute temperature; ©, 
potential temperature; m = &/cy, ratio of 
specific heats; 
a =m q/p, coefficient of piezotropy for 
adiabatic changes 
s log (p™/q) = log © + const, a linear 
function of specific entropy 
= log [(1/p") Vq/r] = 3 [log (t/rp”) — 5]; 
s’=s' + log Vs 
No = Vgs:, coefficient of hydrostatic stabil- 
ity 
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f) k; N; P=2n/N, wave-number, orbital 
angular-frequency and orbital period of a 
sine wave in x-direction, respectively 


g) Dx, D'y, D'z, D'Y, D'Z = Vqr (Dx, 
Du. De DYNDZ) 
ap = Vrla ap; Op = O'p| Vs 

h) 9, latitude; tan ¢ = s,/s., slope’ of the isen- 
tropic surface 
D — + «, angle of Y with equatorial plane 


ı/R = o,, meridian curvature of the level 
surface 

w, absolute angular velocity about earth’s 
axis 

w, vector of magnitude © parallel to earth’s 


axis 
@», earth’s angular velocit 
0» s yi 

u =r (© —w;), eastward velocity rela- 
tive to earth 

f  =2 @, sin y, Coriolis parameter 

F =2wsing, vertical component of 2 W 

F, =2o sin ®; F, = 2 w cos ®, compo- 
nents of 2 W along Z and Y, respect- 
ively 

1 i =— (pz/q) SEE Fir Oz 


€ = F3/(N2? —N®); F’ =F, +r wy 

F,’ = F,—r wy, component of absolute 
vorticity normal to the isentropic 
surface 

G, = FF, — N?;@ =G,+ & (Fr wy — 
— N?) 

EN Ne 

B = (1/G) [FR (F—uy)y + FyN?| + (F/s,) 
(aFu; + s;y) 


4. Specification of the motions 


The undisturbed motion is assumed to be a 
zonal flow of limpid and unviscid air, defined 
by continuous meridian distributions of © and 
of two of the parameters p, q, T, ©, a, s defined 
in (3.e), these three distributions being related 
together by the motion equations, given in 
next section. 

The perturbation is simply assumed to be 
adiabatic and laminar. 


5. Vectorial perturbation equations 


If use is made of the cylindrical coordinates 
L, r, h defined in (3.d), the mechanical 


equations for an unviscid fluid motion may be 
written as follows: 


Bi fie DEN Mie 
2 ODE 
al 2) Hip 


D?r DEN? ” 

Dr? (Zi) f+ B= 0 | (5.1) 
D?h 

prot =O 


Applied to the undisturbed motion, they give 


Pr = 93 (p,lg) + E, = row; (p,/q) + Ex =0 


If we differentiate the last two of these relations 
with respect to h and r respectively, and then 
subtract, we get 


(1/4?) (pn — up.) =2@r@, (5.2) 


or, introducing the parameter s defined in (3.e), 


(5-3) 


and these relations are two equivalent forms 
of the so-called geostrophic equilibrium con- 
dition. It follows from (5.2) that the un- 
disturbed motion is barotropic if, and only if, 
@ is independent of h. 

Instead of (5.3) we may write, in vectorial 
form, 


Vs x (Vp/q) = @w-rVo)x 


The mechanical perturbation equations are 
readily deduced from (5.1) by applying the 
operator 9 to the three equations of this 
system. We thus get, using (3.3), 


(1/9) (pr — Spr) = 2@ ror 


(5.4) 


aD (CDE att) os ate | 
I(r or) tear? ö 

Dr DL\? 19 ae 

? Dr (2) DT 4 
Dh 12 24 

a ain Pr a ] 


But in the undisturbed motion we have 


DS DE a ah Dire Dh DL aan 
DA Dee LANDE ADI 
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whence, using (3.2) and (3.3), and since 
rDL = Dx, 


D2L 
a ( 2x) =D) 2 (« a) =r2D—— + 


DENG IDE) DEN? DE Dr 
Dr? D? D 
ee. Sts DE eh ae 
+ oD Ppp DE t+ 2er5 Dr 
D2 
= (25 De | 205 Dr); 

D?# D?r D? D*h D 
Ina Ppa pal’? pas pal” 
DL\?2 IDL 
(5) de du 

IDE. ER D at 
= EE polo Die) 
201(D Di Do) 20(5Ds r 2) 
We thus arrive at the system 
D? D rue | 
De re nT ee ER 
D? D a3 Ia 
ana By pee rh Ur: I 
pel" 20(5. 5e Da) += 5 ap 
f 
LO 
| 
D? ey) Pn 24 
ee am 
pe # : a q 4 2 ] 


which is equivalent to the vectorial equation 


D? 


pe PM + 2w X > DM — (rDo) x | a 
0 
yo 7 0e, 45.5) 


Let us now consider the continuity equation 
and the thermal equation. For the motion 
itself they may be written respectively, if 
the air is treated as a perfect gas, 


Bee 
Dt 


1 Dq 


. Ds zu. 
Do x 


di aa DE 
ah SB 


and the corresponding perturbation equations 
may be written at once from analogy, if a 
proper choice is made of the correspondence 


between the fluid elements in the undisturbed 
and in the disturbed motion, as follows: 


(5.6) 


From the last one of these equations we get, 
according to the definition of s, 


D 
RER ANSE A PE: Mme 
q q P 


ip 


and if we substitute in (5.5) and the first 
equation (5.6), we get the desired system of 
vectorial perturbation equations as follows: 
D? 


E DM — (rDo) x| 
{ 
f 


Mb ages (V m TE) ap = 0 
q 


div DM +m—2+% dp =o 


6. Perturbation equations referred to 
isentropic coordinates 


Once a coordinate system is chosen, it is 
easy to derive from (5.7) a set of scalar equa- 
tions in which the unknowns are dp and Dx, 
Dx,, Dx, projections of DM on the directions 
of the coordinate lines: these equations are 
obtained by substituting to (5.7.1) its pro- 
jections on these three directions, expressing 
Do, Ds and Dp in terms of Dx,, Dx, and 
Dx, by means of (3.1), and likewise for 
div DM by means of Gauss’ divergence theo- 
rem for instance. We thus get a system of 
four differential equations which are obviously 
linear and homogeneous with respect to the 
four unknowns, and therefore form a com- 
plete system of perturbation equations. 

In the present case of an adiabatic motion, 
it is found that the simplest equations are 
obtained if reference is made, not to the 
usual coordinate lines formed by parallels, 
meridian level-lines and verticals, i.e. tangent 
at each point to the x, y, z-axes, but to the 
similar system tangent to the x, Y, Z-axes 
(see (3.d) for definitions). A coordinate system 
corresponding to such coordinate lines may 
be conveniently called an isentropic system. 


ADIABATIC PERTURBATION EQUATIONS FOR A ZONAL CURRENT AI 


If isentropic coordinates are used, we have 
MODS = sz.DZ..since sy = 0, and the 
system of perturbation equations equivalent 
to (5.7) may be written as follows: 


x) 
m Gi 0, 
p P 


? Dt 
2 
+(e +N})DZ "(2 me?) apo, 


Dr IZ 


0 : à Pr | ry 
= Dx + (4 en ®,) DY + 


| 
{ 
D 
ren Bora DY:+ 
ad LZ. ite CE 
+(2- Lex a + ®y) DZ+ “apo. 
) 


(6.1 


(see (3.h) for symbols). 

The quantities ®z and ®y which appear in 
the last equation are equal to the curvatures 
of the coordinate lines tangent to Y and Z 
respectively, and the following relations can 


be established: 


Dz = as gy = nn © _8¥ sin ¢ (6.2) 
(Oy = py + ey) 


On the other hand, the geostrophic equili- 
brium condition (5.4) takes the form 


er 


—*—sz =Frwz + Rroy (6.3) 


whence we conclude that the coefficient of 
DZ in the second equation (6.1) is equal to 
the coefficient of DY in the third one, and a 
glance at the other coefficients of the system 


reveals that most of their terms are either 
symmetrical or anti-symmetrical with res, ect 
to the principal diagonal. (The same result 
holds with any other meridian coordinates.) 
The symmetry is however incomplete be- 
cause of the presence of the parameters 1/q, 
ry/r, rz/r, Dz and Dy, but it is easy to get rid 
of the first three of them, by using as new 
unknowns the quantities D’x, D’Y, D’Z, Op 
defined in (3.g). We thus get the following 
system, equivalent to (6.1) and obviously 
simpler: 


Ue gee DAR, LA 
RAD N ess alge Owe 


D 
— hm Det Broy DY + 


D? 2 ’ 0 ANS oe 
+ (+ Nt )Dz jr (2 +2) 


Ra RTE. 
EE (£ 


5 | ,) DZ adp-= 0. | 


(6.4) 


sy de) IDE + 


where s’ is the parameter defined in (Zc). 

It is even possible to go further, and to get 
rid also of the term ®z by taking as unknowns 
the quantities 


D’x/Vsz; D'Y¥/Vsz; Vsz D'Z; O'p|V/sz 


The symmetry of the corresponding system 
of equations is practically perfect, since Oy is 
always small compared to sz. However, in 
view of the following developments, it is 
more convenient to retain the system (6.4) 
itself, which will be therefore our fundamental 
system of perturbation equations. 
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7. Zonal waves and wave equation 


A system of differential equations such as 
(6.4), the coefficients of which are independent 
of L and t, is satisfied by solutions of the form 


(DO DYI D's 97) = (Dx, D'Y, DZ, 
Ip) exp i (nL — n't) 


where n is a real positive integer, n’ a real or 
complex constant and D’x,, D’Y,, D'Zy po 
are functions of the meridian coordinates 
only. The corresponding perturbations may 
be conveniently referred to as zonal waves, 
for which we may speak of a wave-number k 
and an orbital angular-frequency N, defined by 


(7-1) 


so that both of these parameters are, as a rule, 
functions of the meridian coordinates. 

In order to get the perturbation equations 
for such waves, we only have to make in 
(6.4) the substitutions 


k = nfr;, N=no—n 


die ik; DIDI=IN (72) 
If this is done, we observe that in the first 
three equations no differential operator figures 
in the coefficients of D’x, D'Y and D’Z. We 
can therefore easily solve these three equa- 
tions for D’x, D'Y and D’Z, and substituting 
the results in the last equation we get a differ- 
ential equation for 0’p only, which we shall 
refer to as the wave equation of the perturba- 
tions. 

In the wave equation, the differential opera- 
tors 2/0x and D/Dt may be reintroduced in 
place of ik and iN respectively if these para- 
meters have been first removed from the 
denominator, and the differential equation 
obtained in this fashion may be more con- 
venient than the wave equation itself to give 
the solution of a perturbation problem when 
the boundary conditions are not simple. 
However in many cases the solution is more 
readily obtained by using the waves as ele- 
mentary solutions, according to Fourier’s 
general method. 


8. Simplification of the wave equation 


The wave equation generally appears with 
a rather complex form, but it happens that 
it can always be greatly simplified if account 


is taken of the order of magnitude of the 
parameters. It is important to note here that 
the simplification must be made only once the wave 
equation is completely derived, and not in the 
course of this derivation. We shall see at the 
end of the paper to what kind of errors one 
may be led by the non-observation of this 
fundamental rule. 

We must first observe that, as a rule, the 
scale of space-distribution of d’p and of all 
its partial derivatives is smaller than, or at 
most comparable to, that of any parameter Q 
involved in the wave equation, which we 
shall refer to as a fundamental parameter. Other- 
wise the portion of atmosphere under con- 
sideration would appear as a surface of dis- 
continuity for Q, as far as the perturbation is 
concerned, and this would not agree with our 
assumption that the motion is everywhere 
continuous. 

If P denotes the variable 0’p or any of its 
partial derivatives with respect to Y or Z, 
we may thus write, as a rule, 


d 2: ) 
QYPE QP; OrPZ Q= Poe 


OZ 
whence 
Ya) 0 d a) 


In other words, the order of magnitude of a 
term of the wave equation is not changed by an 
arbitrary permutation between the parameters and 
the differential operators. 

The relations (8.1) and (8.2) still hold, as a 
rule, if Q is replaced with a function Q’ of 
fundamental parameters, provided that Q’ is 
not simply the difference Q,--Q, between two 
approximately equal quantities (Q, + Q,), or 
a function of such a difference. Therefore if a 
fundamental parameter may be put in the 
form 


Q = +. Q, with Qt <Q; (or OR, 
we can write 


dy Oe 9 


d 
—3y QoP 


and likewise with the operator 9/9Z. We thus 
get 
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d 0 0 0 
ay VP Ht Fy Qo Ps 55 QP H 5G QP (8.3) 


a formula which is very convenient to simplify 
the coefficients of the wave equation. 
Considering now two arbitrary funda- 
mental parameters Q, and Q,, we can easily 
verify, by assuming successively Q, > Q, 
and, Orr that, as a rule, 
Q, + Q~ Q,— Q; (8.4) 
Q, Q = Q == Q > (Qı es Q,)? (8.5) 
and these relations still hold if Q, and Q, are 
replaced with arbitrary symbolic functions 
of fundamental parameters and differential 


operators applied to the variable 0’p, so that 
we have, for instance, 


d oO 
(4 ‘i Q) PE. (3 Qi) p, 
DO Pa: Lo: |? 
1 ANT = 1 2 )Y2 ? 
2? 


0? 0% 
Ren 2 are 
Q, Q, aVaz? = (ei IY? = 03 7) P, ete: 


By these formulae the different terms of the 
wave equation can be compared. In particular, 
it follows from (8.5) that 


1 Q,Q, < Q? + Q3, provided À < 1 (8.6) 


another important formula which will be 
largely used to simplify the wave equation. 


9. General approximate wave equation 


If we solve the first three equations (6.4) 
for D’x, D'Y and D’Z we get, in the case of 
a wave, the following symmetrical expressions: 


oe 1 Fr wy : | 
De=z| (+ a) 


TOME 
cl le +) | 
Pee ] ! 
a 7 
ee | “) Ex 
He EN F, F,’ 
Dale Line N?—N? 


(see (3.h) for symbols). 


If we substitute from (9.1) in the last equa- 
tion (6.4) and observe that, according to 
(7.1), we have 


ky/k =—ry/r; Ny = kroy 


(9.2) 


we obtain, for the wave equation, 


(41 + A, + A3 + Aj) d'p =0 


(9.3) 


= a 


| Gi 0 
G'(N? — N3 \aZ 


3 Fk [Gy ur 
GIN: Gi JER r Sz 


F,Nk ig +) zu 


1 
dy). 
F, F,' ig AS 

"Gm: — N5 (on a s) 


Let us note here the important result that the 
term in k?/N? is exactly cancelled out in the 
expression for A, in virtue of (9.2), at the 
same time as the terms in 9/9Y cancel out each 
other. 

Let us now refer to the current values of 
the parameters, as observed in the meteoro- 
logical atmosphere. We have in cgs units 
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sg~~ 1077; —pz/q Hg + 108; F, #f 1074; | 
Ee HO ae TON: (Oy, aud uy =f a ane 
uz = 107 3. Selon. 1/R~ 10- os Dull aa | 


(9.4-1) 
whence 
Fy Er WS So, To | 
eh Sy Or Sy te $7; IN te NG TO 
PALF ae GH Ct Ch We “| 


(9.4.2) 


Using these relations, we can prove that, in 
the wave equation, 


i) G, and F, may be substituted every- 
where to G’ and F’ respectively, and e’ may 
be omitted in the second term of A,. This is 
obvious in virtue of (8.3 
successively. Q, =G; nd "Qi='G"; -Q, = Fi 
ADOROENES Ox stand GE, 


ii) G,y and F,y may be likewise substituted 
to Gy and Fy in the second expression for Ag. 

We have to prove that we may neglect 
the two quantities 


i MER es 

ay == (ee GN (G G,)y; 
i,k 

ET GN (dr F)s 


From the definitions of G’ and F’ we get, 
since krwy = Ny I = N (d/0Y), 


i, TAF ; Klar ee 
m= Gan (FC —G) ay = & 
EEG NE Sead” N 
Ca re Ayo 
I k do £” on 
ee F’ ~ gh 
Zorn Ay = Gap 


whence, since F’F,’ = G’; 
and according to (8. 5); 


ae ee | LE ee 
ae Naya DY) Gs LOT Spa 


ee A 
di Se ee 
2 G! a y2 < 1 


FUN Sr Gare er, 


ii) Py may be omitted everywhere. 


) where one makes. 


Owing to the relations Oy = py + ey; 
py < sz’; ey Te (A/AY), we have only to 
prove that we may neglect the three quantities 
L , F,Ne ok 
ZI CHND. EE Na 


€ a) a) 
NEN? JXNGZ. 


FAP, € Ja sit 
GN: — NS ay \ay ' *¥ 
Since we have, by (8.5), 
AS = G k G' (+ = s,) at 


it follows from (8.6) that these terms are 
negligible compared to A, provided 


eV = <1; € FN 
NI NE nl NME NE 


€ - 


À 


€ I 
2 + N 
and these three Er ad are actually satisfied. 
iv) The remaining terms of A, and A, are 
negligible compared to A,. Since szy/sz = 
> d/9Y, this is readily proved by putting in 
(8.6), successively, 


I I ) 
Qu. = by Ole ae ae 
1 VG 2 VN? — N° 2Z € GE re 


LRU x); Ne en 
VG (Ni— N94 


for in both cases we have A <= ı 


v) N, may be substituted to N, in the 
expression for A,. Let us put 


ref TE 
LOL, 4 Mens 
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If we refer to the definitions of N, and N, 
and to equations of motion, we find that 


E1 = & + ul — 2), with 


szr cos ® 


: Fy raz 
KT mm’ B 
1 


N? — N? 


whence, by (8.2), 


0 IL, 
SZ SZ d A 
€ — — - 2) 28,0, — 
252 Nic NB E Soyo 
SZ 
—€& + Co € Gy 
SZ 


in- Other words, Ci +-iCp. 


We thus arrive at the following simplified 
wave equation: 


[= d 2, SZY T d 
PR A LG; OY 


: Boe. 2 ry SZY 
B= F(Z tay th 4 
This equation can be transformed in several 
ways. 

We can first prove, by the same methods 
used for above simplifications, that G,, Fy, 
sz, sz and 9/9Z may be replaced everywhere 
with G, F, s., 5 and 0/dz, respectively. That 
means that we may disregard the small de- 
parture of Z-direction from the vertical, as 
we could expect it from the smallness of all the 
horizontal gradients compared to the vertical 
ones, and that furthermore rwy may be re- 
placed with uy. 


On the other hand, if we refer to the de- 
finition of s’ and to relation (6.3) we get, 
since sy = 0, 


= . (Fyroz + Fyrwy) 


and, according to (9.2), 
Gy =F (F — uy) y + Fy (F— uy) — 2 Nkroy 


whence we get, after a few obvious simplifica- 
tions have been made, the following form for 
the wave equation: 


[ (+ 2e) k? d ; S:Y\ I 
© Nemo nn Rules 


ARS: Vu A, I 
IN dz 2 8 N? Im 
ig Oe k | 
(2: x) ra an] PO 
(9.6) 


where B is the quantity defined at the end of 
G-h). 

Another practically equivalent form is 
obtained by introducing the variable 0”p and 
the- parameter }s” [see /(3.f) and (3.9], “as 


follows: 
| A 2 Fuy k? 0 “ I 
de av IG 
ak NG 
A N? N? 


Lf if a) 
z N? \2z 
ta 2] apo 


ie jt 
(2 +H") GN 


[In (9.6) and (9.7) roy may be substituted 
to uy, and the same substitution may be made 
in the expressions for G and B.| 

We observe that the vector 2 W only figures by 
its vertical component in the wave equation 
(9.6) or (9.7), in agreement with the current 
assertion that the horizontal component of 
earth’s rotation may be always neglected in 
atmospheric motions. 


(9.7) 
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Let us now compare the orders of magnitude 
of the different terms of the wave equation 
(9.5). We have, when dealing with operators 
applied to d’p, 


ke 3 . SN Le e +8) > 
(es DATE EEG, CAP) 


I A Ne 
- a led & »)]. 


d I 7 “8 Er 
Brise EN IN ze 
I a LEN 
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This formula will be very useful in the next 
sections. 


10. General approximate formulae for the 
vibration components 


In the case of a wave, the vector DM may 
be conveniently called the vibration, and its 


= 


components are given as functions of d’p by 
formulae (9.1). Let us now show that these 
formulae can be simplified like the wave 
equation, as far as the magnitude and direction 
of the vibration is concerned. For this purpose, we 
must observe that any term of the expression 
for D'x, D'Y or D'Z may be omitted if it is 
small compared to any one of these three 
components. 

i) The last terms of D’x and D’Y are 
negligible compared to the other terms of 
these components. 

This is readily seen by using relation (9.8), 
whence we deduce that 


EN Mas ee F,N 
EEE = d >= en 
N? — N° > z= > P< JGIN EN) 


Fe a ; : 
(24+) (K+ ay 7 sy) ap = 


Ven a) ; 3 
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PAF,’ a AE ble 
NN (a2 : 4) fu VG (N2 —N?) 


jee d ; ! 
(+5 (+ + ss) ap = 


) 
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ii) The first two terms of D’Z are negligible 
compared to the first term of D’x and the 
second term of D’Y, respectively. 

This results from the fact that 


rN FRY _ 
N: = N2 N2 = 2 = 


1 


< ı and 


mic 


ii) 2/0Z and $ may be replaced with 
010Z and sz, respectively, in the last term 
D we. 


Similar proof as above. 
iv) D'y and D’z may be computed by the 
same formulae as D’Y and D'Z, respectively. 
This is obvious from the relations 
D'y = D'Y cose + D'Z sins; 
D'£ =) DZ € DY Bins 


Finally, if we simplify the terms of (9.1) 
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according to (9.4), we arrive at the following 
approximate formulae: 


ee E 1 Froy 
Dx= &[—(:+ am) + 

J'EN iby? ; j 

as x (y+ 6) 2p 


By Yor Dy = 5 [+ (or - )] op 


Iron Di = 


I a 
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I 0 . 
Ni — N° Ë t «) ge 


(10.1) 


+ d'p or — 


We observe that, just as in the wave equa- 
tion, the vector 2 only figures by its vertical 
component. 


II. Classification of waves according to 
orbital period P 


Up to this point, we have made no assump- 
tion about N. We shall now see that further 
simplifications appear if this variable is assumed 
to have a definite order of magnitude com- 
pared to the parameters N, and f. 

Since we have always N,/fS 10%, it is 
obvious that we are always in only one of the 
three following cases: 

i) N= N, > f (P < 2 hours), short-period 
waves 

DEN N and = foi hours < P < 5 
pendulum days), middle-period waves 

ii) N<f< No (P > 5 pendulum days), 
long-period waves 

According to these definitions, the width 
of the middle-period waves interval in the 
frequency spectrum becomes larger and larger 
as latitude decreases, and at the equator there 
are no long-period waves at all. 


i)  Short-period waves 


From the assumption N > f, it can be 
proved, by the same method used in preceding 
sections, that the wave equation reduces to 


E 9 dal 0 ; 
Name À ay moi He elle 


d st) is! d “igs, news| 34 
22 “NN: | * «| Fe 
(11.1) 


On the other hand, the expressions (10.1) 
for the vibration components may be replaced 
with 


i 


/ hate / I 7] 2 
I 9 A 
(2 = s) ap (112) 


N?— N° 
We can see that now neither the vector W nor 
its ascendant figure in the formulae, so that 
we can obtain the same formulae by assuming 
from the beginning that w = o. 


D'z= 


ii) Middle-period waves 


If N<N,„ we can simplify the wave 
equation (9.6) or (9.7) by omitting N? beside 
N?, but a more interesting result concerns the 
vibration components. 

The relation (9.8) gives 


I 0 2 be eee 
N? — Ni (2 +4) EN 


F 0 / , I 

(24+) (k+ B+ %) op<g: 
F 7 ; _ 
(+3) (b+ 3 + 5) 27 


whence we conclude, by (10.1), that 


Daa Dx ior, Dy 


In other words, the vertical component of 
the motion is negligible compared to its 
horizontal components. This is in agreement 
with the observed fact that if the scale of an 
atmospheric motion is such that its vertical 
acceleration is negligible (an assumption 
which corresponds to the condition N < Np) 
this motion is always mainly horizontal. 

Let us observe here that, in low and middle 
latitudes, some of the middle-period waves 
satisfy both of the conditions N < N, and 
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N > f, so that the special properties of the 
short-period waves also belong to these waves. 


iii) Long-period waves 


For these waves, we have of course the 
same fundamental property as for the middle- 
period ones, namely the predominance of the 
horizontal components of the motion, but 
further simplifications result from the condi- 
tion N < f: 

In first place the expressions for G and B 
reduce to 


G=F(F—uy); 
B—F [= uy)y i aFu; + =] 
F— uy Se 


and secondly the horizontal components of 
the vibration may be computed by the ap- 
proximate formulae 


i if ‘Oy, , (0 s)] À 
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12. Wave equation with constant 
coefficients 


a) Let us first consider the case of a short- 
period wave, with 


© — const, whence N = const 
2 
N, = const; 2/dy = 10-8; 9/92 1075] (12.1) 


We thus have sy’ < 9/0y; s,’ < d/dz 


so that the wave equation (11.1) reduces to 
k? I 2 I 0° 4) 
; — — —_-__. __ __q]} '»=0 
N: Mop N®—N? gz? oP 
(12.2) 


and all its coefficients may be considered 
as constants, as well as r and R, since 
their scale of space-distribution is, according 
to (12.1), large compared to that of the 
perturbation. We may therefore take, as 


particular solutions, the pseudo-plane waves 


defined by 


d'p = dp, exp. i (kx + hy + k’z—n't) 

(12.3) 

where x =rL; y = Ry; z is the altitude, 
and d’p,, k!, k’ are three constants. 

From (12.2) it follows that k, k’, k” and 


N must be connected together by the 
relation 


(RAR PEN + KIN NS) ae 
or 
a Ni—(K? +a N2) N° + N2 K?sin?y =o 


where K = Vk? + k2 +k” is the geo- 
metrical wave-number, and y is the slope 
angle of the wave surface (tan y=k”/K). 


a.1) If k’ and k” are real (internal waves 
with uniform amplitude), it results from 
(12.1) that KY 1075, whence a N < K? 
and VaN, sin y < K, so that the above 
relation may be written as well 


(a N?— K?) (N?— N32 sin? y) =o 


We thus get two types of waves: 
i) If aN?— K? =o, the orbital phase- 
velocity U = | N/k| has the uniform value 
U = 1/Va, the sound velocity, and since 
N > N, the relations (11.2) give us 

Dx Dy Dz 


Ee, Oke LS 


a proportion which means that the vibra- 
tion is longitudinal. The corresponding 
waves are obviously nothing but the clas- 
sical elastic waves. 


ii) If N®— N? sin? y =o, whence | N|= 
= N, sin y, we get from (11.2) 


kDx + k’Dy + k"Dz =0 


in other words the vibration is transversal. 
We have proposed the name internal gra- 
vity-waves for the corresponding waves. 

a.2) If k’’ is a pure imaginary number, 
k’ being real, the waves are no longer 
internal waves with uniform amplitude, 
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they are surface waves with an amplitude 
which decreases exponentially upwards or 
downwards. For instance if we assume 
N =" No, and aN? = 1071%,. we get 
k''2—= — (k2 + k’?), and the corresponding 
waves are of the type of the classical 
gravitational waves observed near the sur- 
face of deep water for instance. 


Let us then consider the case of a long- 
period wave, with 


@ = const, whence again N = const | 
A} 9 

N= canstd kiss 207%: d/dy ~ 10°] (12.4) 

Since = const and s, + const, we have 


G=F; B+ Fy +f, 


and the wave equation (9.6) reduces to 


5 


If the study of the motion is restricted 
to a narrow middle-latitude zone of say 
20° or 30° in width, all the coefficients of 
this equation may be taken as constants, 
as well as r and R. We may therefore 
again consider particular solutions of the 
form (12.3), with now k’=o, and 


0? k 


RL) d'p FLO 


F2? 
N 


dz? 


k2 + (F2/N2) k’? — (k/N) fy =0 


If we consider the orbital phase-velocity 


U, =—(N/k) in the x-direction, we get 
from this relation 
U, =—fy/[k? + (F/N2)k'? (124) 


and if, as a result of proper boundary 
conditions, the term in k” is small com- 
pared to the term in k, this dispersion 
equation reduces to RossBy’s formula for 
long atmospheric waves. 

If in (b) we replace the condition w = const 
with the two following ones, which are 
practically compatible: 


absolute vertical vorticity = F — u, = const; 
bye 6 


we have F— uy = const, whence B = 0, 
and the wave equation becomes 


Frey : Sy =| ’ 
RE 
| Br ( eg. 


4—906620 


13: 


2 
en (3 


49 


If the study of the motion is again restricted 
to a narrow zone, the coefficients of this 
equation may be taken as constants, and 
we may again consider solutions of the 
form (12.3), with k = o and 


Li Et N 
ay F 


F — uy 
Considerations about some incorrect 
results recently published 


L'"’ ak 


k 


Some authors have considered that, since the 
large-scale motions are always mainly hori- 
zontal, it is permissible, in the case of 
long-period waves, to disregard the ver- 
tical motion in the perturbation equations, 
a procedure which obviously facilitates the 
derivation of the wave equation. It may 
be useful to show how far the wave equa- 
tion obtained in this fashion is incorrect. 

If we neglect the vertical motion, we are 
led to replace the perturbation equations 
(6.4) with the following approximate 
ones, referred to x, y, z-axes: 

D? D 
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whence, by eliminating D’x and D’y, we 
easily get the wave equation 


s,) D'y + ad'p=o. 
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Let us now consider the special case where 


Wy = 0, whence N, = 0, and 
N? < F®; k = 10-8; d/dy = 10” 


With these assumptions, (13.1) reduces to 
the following dispersion equation: 


(: Fi) k2 zn =) 


with now 
By + fy ki f (afu +, sy) 


This is a second-degree equation in N, the 
roots of which are imaginary if 

Be x 
4 | Jus, | 
so that the above derivation leads to un- 
stable waves in a westerly current if the 
potential temperature decreases towards the 
pole. It is clear that this result, which has 
been given as an explanation for the origin 
of some of the extra-tropical atmospheric 
waves, is due to the presence, in wave 
equation (13.1), of the term in k?/N?, 
which has not been cancelled out in the 
derivation of this equation, precisely be- 
cause the vertical component of the mo- 
tion has been neglected. Since such a term 
does not exist in the correct wave equation, 
as already emphasized, we must conclude 
that the existence of unstable waves in 
the westerlies cannot be justified in this 
simple fashion. 

By this example, we also understand to 
what kind of errors may lead the applica- 
tion of the method under discussion. 


as, = ON and Re 


Another simplification which is frequently 
made is that, according to what is assumed 
in the so-called parcel method in the 
theory of hydrostatic stability, the local 
change of pressure dp is negligible, so that 
one may put dp =o (or d’p =o) in the 
wave equations. 


If we do that in equations (6.4), and if we 
only retain the first three of these equations, 
we get, by eliminating the unknowns 
D'x, D'Y, D’Z, the following condition 
in the case of a wave: 


G +f F(F—rwy) — N?®=o (13.2) 


which determines the value of N at any 
point. 

The corresponding waves, in the assump- 
tion that they can exist, are nothing but 
the so-called isobaric waves, so often consi- 
dered in theoretical meteorology, and the 
condition (13.2) shows that these waves 
would be stable or unstable, according as 
roy is smaller or larger than F. If for in- 
stance w = const, we get N = + F, and 
if © = w, the orbital frequency must be 
equal to the Coriolis parameter f. 

But it is easy to see that such waves cannot 
exist except for particular undisturbed 
motions. For N is already connected to w 
by the relation N = nw —n’, which is 
compatible with (13.2) only for a special 
space-distribution of w. Furthermore, what- 
ever may be this special distribution, the 
isobaric waves are necessarily stable, for 
nw —-n’ can bea pure imaginary everywhere 
only if © is a constant, and we have already 
seen that in this case N is real. 

We can therefore conclude that the con- 
sideration of isobaric waves is only of poor 
interest in the case of a distorted current and 
cannot explain the apparition of unstable 
waves. 


14. Conclusion 


This paper was written with the intent, not to 
give the solution of a particular perturbation- 
problem related to a zonal air-current, but 
only to secure the first step in the way to such 
a solution, namely the derivation of the wave 
equation. The perturbation problem is thus 
reduced to a purely mathematical one: the 
resolution of the wave equation under proper 
boundary conditions. 

This resolution is certainly not an easy task 
for most of practical problems, but it is clear 
that it will be always possible to manage it, 
by using modern computers if necessary, and 
we have proved in the preceding section that, 
in fact, most of the disagreements between the 
results obtained by various authors are to be 
imputed to incorrect derivations of the wave 
equation. 

We think it therefore not superfluous, in 
view of generalizations to non-zonal or non- 
adiabatic motions for instance, to emphasize 
here that, when deriving a differential equa- 


ADIABATIC PERTURBATION EQUATIONS FOR A ZONAL CURRENT SI 


tion, the simplifications resulting from the order 
of magnitude of the parameters must be made only 
at the very end of the derivation, and not in the 
course of the computations. According to this 
principle, the approximate meteorological 
laws generally applied to large-scale motions 
(vanishing of vertical acceleration or velocity, 


vorticity theorem or conservation of vertical 
vorticity, etc.), though very useful to give the 
distributions of some of the hydrodynamical 
quantities when some others are known, 
must never be used as starting relations in a 
sequence of computations where differentia- 
tions are involved. 
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Sur la circulation transversale associee A un courant 


atmospherique 


Par JACQUES VAN MIEGHEM, Université de Bruxelles 


(Manuscript received 14 December 1949) 


Abstract 


A criticism of the interpretation of the circulation theorem is followed by an investigation 
of the different factors which contribute to a transversal circulation in a horizontal current. 


1. L'étude de la circulation transversale 
associée au et-stream» [5] a conduit aux 
notions de circulation directe et de circula- 
tion inverse, au sujet desquelles il faut, croyons- 
nous, faire certaines réserves. Une circulation 
est dite directe lorsque son sens est celui pre- 
scrit par les solénoïdes isobares-isostéres; elle 
est dite inverse dans le cas contraire. On admet 
qu une circulation directe résulte des solénoïdes 
et ne se maintient qu'aux dépens de ceux-ci, 
tandis qu’une circulation inverse est supposée 
engendrer des solénoïdes isobares-isostères. De 
plus, confondant les notions «sens de circula- 
tion» et «sens du cycle thermodynamique dé- 
crit par une particule d’airm, on affirme qu’une 
circulation directe transforme de l'énergie po- 
tentielle en énergie cinétique, une circulation 
inverse, de l'énergie cinétique en énergie po- 
tentielle [3]. 

L'interprétation du théorème de la circula- 
tion de V. Byerknes à la lumière d’un principe 
de causalité ne peut conduire qu’à des con- 
tradictions. En voici un exemple récent. D’après 
P. RAETHJEN [4], le gradient isobare y, de la 
température crée une circulation qui rétablit 
le gradient défini par le taux d’accroissement 
du vent suivant la verticale. Lorsque le gradient 
Yp est trop fort (c'est-à-dire, supérieur au gra- 
dient isobare 1°, défini par le taux d’accrois- 
sement du vent réel suivant la verticale), l'air 
chaud s'élève et lair froid descend (circulation 


directe), les variations adiabatiques de la tem- 
pérature de l'air qui en résultent, entraînent 
une diminution du gradient y, tant que la 
valeur limite J’, n’est pas atteinte (condition 
d'équilibre de Marcutes). Inversément, lors- 
que le gradient y, est trop faible (y,< I>), 
lair chaud descend et l’air froid s’eleve (circu- 
lation inverse) jusqu'à ce que le gradient y, 
atteigne la valeur limite 7°. E. PALMEN (3) 
par contre, prétend que les forts gradients de 
température que l’on rencontre dans la surface 
de soo mb où le «jet stream» contourne un 
creux, sont dus à une circulation inverse 
susceptible de créer des gradients de cette im- 
portance, laissant sous-entendre que les faibles 
gradients thermiques dans les dorsales de la 
même surface isobare, doivent être attribués à 
une circulation directe. 

Ainsi, pour P. RAETHJEN, les solénoïdes 
constituent l'agent moteur de la circulation, 
tandis que pour E. PALMÉN, la circulation est 
créatrice de solénoïdes. De telles contradictions 
sont inévitables dès que l’on fait appel à un 
principe de causalité. Qu'il me suffise de 
rappeler que les équations de la mécanique 
expriment des relations entre certaines gran- 
deurs, mais qu’elles n’établissent aucun lien 
causal entre celles-ci. 

2. Le théorème classique de V. BJERKNES 
exprime que dans un fluide barocline, la cir- 
culation absolue est une propriété non con- 
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servative du mouvement du fluide et cela par 
une relation entre le champ du mouvement 
et le champ des solénoïdes isobares-isostères et 
non pas par un lien causal entre ces deux 
champs. L’accélération de la circulation ab- 
solue le long d’une courbe fluide et le nombre 
des solénoïdes isobares-isostères, que cette 
courbe entoure, varient simultanément dans le 
même sens, sans que l'on puisse affirmer que 
les solénoïdes sont la cause ou la conséquence 
de la circulation. Enfin, le sens de cette circu- 
lation est univoquement déterminé (sens so- 
lénoïdal). 

Dans un système de référence fixe par rap- 
port à la terre, l'accélération de la circulation 
relative le long d’une courbe fermée fluide est 
égale à la différence de deux termes: le premier 
n'est autre que le nombre de solénoïdes iso- 
bares-isostères embrassés par cette courbe 
(c'est-à-dire l'accélération de la circulation 
absolue le long de cette courbe), le second 
dépend de la rotation de la terre et de la dé- 
formation du milieu fluide à l’intérieur de la 


courbe envisagée. Le sens de la circulation’ 


relative est déterminé par le signe de cette 
différence. Cette circulation peut donc avoir 
le sens solénoidal ou le sens opposé, sans qu'il 
soit pour cela possible d’affirmer que la circu- 
lation relative résulte de l’effet differentiel des 
deux termes susmentionnés, ou inversément, 
que la circulation engendre cet effet par une 
création ou une destruction appropriée de 
solenoides isobares-isostères. 

3. Pour étudier la circulation transversale 
associée à un courant horizontal en un point 
quelconque O de ce courant, il est commode 
d'adopter un trièdre trirectangle dextrogyre 
Oxyz, fixe par rapport à la terre, dont l’axe 
Ox coïncide avec la direction du courant et 
l’axe Oz, avec la verticale ascendante de O. 
Dans ce cas, l’axe Oy est orienté à gauche du 
courant et par conséquent, dans la troposphère 
moyenne et supérieure de l'hémisphère Nord, 
cet axe est dirigé, à la fois, vers les basses 
pressions et les basses températures. 

Dérivons par rapport à y l'équation du 
mouvement suivant Oz et par rapport à z 
l'équation du mouvement suivant Oy; en 
soustrayant les relations ainsi obtenues, on 
trouve 


dv\ 9 (dv; 2) 7) il 
ee . ~ ay Zz) 2e \dt 


du Op dx op MW». 


2) iy dz 


di Oz dy 


Oz 


Comme d’habitude, p et « représentent la 
pression atmosphérique et le volume spécifique 
de lair, vy, vy, v, les composantes coordon- 
nées de la vitesse © de l'air par rapport à la 
terre (on: a3, v1 > 10..et et. O,.0, a) 0 Chats 
& @z les composantes coordonnées de la 
rotation de la terre (on a: &w, = cos y sin B, 
©, =o cosy cos P, ©, =o sing, où y désigne 
la latitude et 6 l'angle du courant avec la 
direction W — E). 
Ensuite, remarquons que 


0 . 
en = — — sin @ cos Bsin B et 
di 1 
on = — = sin y cos:p; 
d’où, 
2 ac ) a = ze sin œ cos B 
1 dy zei Aid) @ ? 
(v, cos B — v, sin B) = — EH COS D. (2) 


u étant la composante zonale du courant % 
(en O, on a: u = vz cos ß) et a, le rayon moyen 
de la terre. à 

Le nombre N, de solénoïdes isobares-isos- 
tères, par unité d’aire, en O, dans le plan 
transversal de ce point, est donné par 


Route du op sg (OT cs 
Ta ee o Toyah 


MW, 


rn) 


compte tenu de l’équation des gaz parfaits, de 
l'équation de la statique suivant la verticale z 
et de l’équation d’équilibre de Margules; dans 


(3), (5 ) représente l’ascendant isobare — y, 
p 


de la température absolue T de l’air suivant 
la direction transversale y, g, la pesanteur et 
v, la vitesse géostrophique de lair en ©. 
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Enfin, substituons dans (1) les relations (2), 
(3) et l’&quation de continuité, on obtient 


dv ™ Wr OW, We 
a Ga) ee Oe a wae Og 
20, de | Ws, cos ß 
Die ray 20 Hrn. (4) 
4 WV, f 
au point O. Rappelons que 7, tepresente le vent 


thermique en ce point. 

L’acceleration de la circulation (ou ce qui 
revient au méme, la circulation de l’accelera- 
tion) le long de la frontiére d’une unité de 
surface dans le plan transversal yz est égale à 


la composante longitudinale rots( de la 
dv 


rotationnelle de l’accélération — de l'air, au 


dt 

oint de convergence de la surface considérée. 
La formule (4) exprime donc l’accélération de 
la circulation, par unité de surface, en un point 
quelconque O, dans le plan transversal de ce 
point; elle s'applique aussi bien 4 un courant 


> 


a 


a composante W (« ag Gi 4 


qu’à un courant à composante 


< pa 22). 


2 

Comme dans la troposphere moyenne et 
supérieure de l'hémisphère: Nord, Tair froid 
est à gauche ct l’air chaud à droite du courant 
v, il se produira une circulation transversale 


EE. ai & 
relative directe lorsque rot, ay lee? (sens so- 


lenoidal) et une circulation inverse lorsque 


rot, FP 0. 


Examinons maintenant la contribution de 
chacun des cing termes du second membre de 
(4) au sens de la circulation transversale relative. 

1° Lorsqu’en un point quelconque O, le 

, : I; 
taux d’accroissement —* du courant v, en ce 
Oz 
point, suivant la verticale z (vertical wind 
shear), est inférieur (supérieur) au vent thermi- 


OV, : Re DA 
que =, en ce point, le terme 2 w, | —*— — 
2 


contribue à une circulation directe (inverse) au 
même point. Ainsi, un accroissement du vent 
thermique (ou ce qui revient au même, un 
accroissement du nombre de solénoïdes isoba- 
res-isostéres) est associé à une circulation directe, 
une diminution, à une circulation inverse [2], 
le gradient vertical du vent réel étant supposé 
constant. 


2° Lorsque le courant v, au point O augmente 
(diminue) dans le sens x du courant, le terme 


Vee : x : B 
20% a contribue à une circulation transver- 


sale directe (inverse) lorsque le courant v, est 
de NW ou de NE (wx < 0) et à une circulation 
inverse (directe) lorsque le courant v, est de 
SW ou de SE (w, > 0). 

On sait que Vintensité du «et stream» des 
régions tempérées varie le long du courant, 
les maxima d’intensité se trouvant un peu en 
amont des thalwegs des creux bariques [5]. 
Par conséquent, entre une dorsale et le creux 
suivant (à VE) où le courant de NW est 


laa , | > ’ 
accéléré, ainsi qu'entre un creux et la dorsale 


suivante (à l'E), où le courant de SW est 
retardé, il se produit une circulation transver- 
sale directe en relation avec la répartition de 
l’accélération longitudinale du courant. 


3° Lorsque dans le courant v, l’air descend, 


da dx 
on a — <o [| monte, — > o) etle terme 
dt dt 


2 Ws da à x à x 5 
— contribue à une circulation directe 


HUE 
(inverse) lorsque le courant v, est de NW ou 
de NE ct à une circulation transversale inverse 
(directe) lorsque le courant v, est de SW ou 
de SE. 

D’après R. G. FLEAGLE [1], lair a une ten- 
dance à l’affaissement dans un jew de NW, 
tandis qu'il a une tendance à l'ascension dans 
un get de SW. Cette circonstance contribue 
À une circulation transversale directe dans les 
deux cas. 


4° Lorsque le coùrant v, a une composante 


7 a 
zonale d’W C >0,—- < P< -eta, >), 
2 2 
dx . À : : 
le terme 2, > contribue 4 une circulation 


; ; ov : 
directe à gauche où (> < 0, tourbillon de 


glissement cyclonique dans le plan horizontal 
xy) et à une circulation inverse à droite (où 
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x 

dy = 0, 

nique dans le plan horizontal xy) du courant v,. 
Rappelons ici que, d’après C. G. Rosssy [5], 

la circulation transversale est inverse au S d’un 


3 3 Ov 
get stream» d’W, où u > o et re 
Y 


Lorsque le courant v, a une composante 


: TT 
zonale d’E (« er NE 


5 
a 


tourbillon de glissement anticyclo- 


37 
ap <= J ety <0), 


Ov, 


le terme 2, —. contribue à une circulation 
2 


à ... fave x 
transversale directe à droite ( 5 Se ) et à 


. Pe WV. 
une circulation inverse à gauche ar = ) du 
oy 


courant vx. 


Il en résulte qu'au S d’un «et stream» d’E, 
il y a aussi une tendance au développement 
dune circulation inverse. 


5° Comme le produit u cos ß est toujours 
positif, quelle que soit l’orientation du courant 


vy, le dernier terme, dû à un effet de latitude, 
contribue toujours à une circulation directe. 
Nous avons décomposé l'expression de 


dv 
or, en cing termes dont le premier et le 


quatrieme sont vraisemblablement les plus 
importants et dont le dernier est sans doute 
négligeable. Dans l’état actuel de nos connais- 
sances, il est difficile de préciser davantage 
l’ordre de grandeur de ces termes. Une étude 
synoptique de la difference entre l’ascendant 
vertical du vent (vertical wind shear) et le vent 
thermique dans le «et stream» nous parait 
indispensable si l’on veut arriver à une connais- 
sance meilleure de la circulation transversale 
associée aux ondulations d’un courant zonal. 
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On the Vertical Distribution of Wind and Temperature over 
Indo-Pakistan along the Meridian 76° E in Winter 
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Abstract 


A mean cross-section over Indo-Pakistan along the meridian 76° E for January and 
February 1946 has been computed. Three jet stream centers are obtained along this meri- 
dian: one in Siberia, one just south of the plateau, and one in the tropics. The effect of 
the plateau gives rise to Siberian and Himalayan jets. A numerical comparison of the com- 
puted mean profile along 76° E with that of Hess (1948) along 80° W over North America 


is presented. 


Computed wind velocity, wind shear and absolute vorticity patterns are discussed, par- 
ticularly with regard to strong anticyclonic shear and very low absolute vorticity found on 
the south side of the Himalayan jet. The spatial distributions of the isentropic field, the 
gravitational stability field, the baroclinity field and the rainfall during the winter months are 
discussed in regard to the 'sense of the meridional circulation around the Himalayan and the 


Equatorial jets. 


Introduction 


Mean vertical cross sections of the average 
zonal wind distribution aloft across the north- 
ern hemisphere from pole to equator were first 
prepared by V. BJERKNES and COLLABORATORS 
(1933). Although at that time the available 
observations were sparse and not homogeneous 
in space or time, the authors’ sections (fig. r) 
nevertheless permit us to deduce that the 
westerlies even in the mean are not a broad 
uniform current. Their momentum tends to be 
concentrated in a narrow meandering river of 
fast moving air—frequently called »jet stream» 
(RossBy and COLLABORATORS : 1947). More 
recently WILLETT (1944) and Hess (1948) have 
offered similar profiles for North America. 
On their sections we find the center of the jet 
stream near latitude 35° N at a height of about 
12 km in winter. We en counter it in sum- 

1 This research was carried on when the author was 


studying at the University of Chicago under the sponsor- 
ship of the Government of Pakistan. 


mer near latitude 55° N, 
intensity. 

Cross sections made for individual days by 
RIEHL (1948), PALMÉN (1948), PALMÉN and 
NAGLER (1948) and mean cross sections pre- 
pared for a certain short period by PALMÉN 
and NEWTON (1948) have demonstrated the 
same feature of concentration of the westerlies 
aloft. As should be expected, the intensity of 
such instantaneous jets exceeds that of the mean 
sections. They also show that frequently two 
maxima are present in the belt of the wester- 
lies, a fact also elaborated on by CRESSMAN 
(1949). A low latitude west-wind maximum, 
separated from the main center of the wester- 
lies at higher latitude, is a prominent feature of 
the tropical circulation over the Atlantic and 
Pacific oceans (1945). 

Until now three different explanations have 
been offered for the existence of the meander- 
ing river of fast moving air. RossBy and CoLLA- 
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Fig. 1. Mean vertical distribution of zonal component 

of geostrophic wind in northern hemisphere for February, 

reproduced from Physikalische Hydrodynamik after 
BJERKNES and COLLABORATORS (1933). 


BORATORS (1947) have advanced a theory based 
on the principle of conservation of absolute 
vorticity in connection with large scale hori- 
zontal mixing. PALMEN and NEWTON maintain 
that the jet stream is associated with the polar 
front while Namias and Crapp (1949) have 
formulated the “Confluence” concept. 

It is well known that the principal polar front 
hardly passes through Indo-Pakistan. There is, 
of course, evidence that disturbances enter the 
sub-continent from the west during the winter 
months. Important low-level concentrations 
of isotherms, however, are rarely noted in con- 
nection with such disturbances. Granted that 
Bjerknes’ cross-section is representative, it is 
difficult to explain the jet stream as based on 
association with the polar front. In view of the 
foregoing discussion it is worthwhile to in- 
vestigate the mean wind distribution over 
Indo-Pakistan in winter with recent aerologi- 
cal data, in order to see whether the more 
numerous observations now available confirm 
Bjerknes’ section. 


Available Observations 


During World War II, a large number of 
radiosonde stations was operated over Indo- 
Pakistan. This network makes it possible to 
construct a cross section along the meridian 
76° E. If Siberian data are added, the section 
extends from Bangalore (13°N) to Omsk 
(55° N). 

From the observational viewpoint, January 
and February 1946 proved to be the most 
feasible period for computation. Data for India 
were available from the published meteoro- 
logical report of the India Meteorological 
Department. Siberian data (soo and 700 mbs 
only) were obtained from the Northern 
Hemisphere Historical Charts of the Air 
Weather Service. The following stations have 
been used. India: (r) Bangalore (12°s9’ N, 
8190 ER), (2)) Poona (18.387N,1 73233..E): 
(3) Allahabad (25°25’N, 81°51’ E) and (4) 
Delhi (28°39’ N, 77°17’ E). Pakistan: (1) Pesh- 
awar (34°2' N, 71°37’ E). U.S.S.R.: (1) Tash- 
kent (41°20’N, 69°18’ E), (2) Alma Ata 
(43°15’ N, 76°55’ E), (3) Novosibirsk (54°48’ 
IN} 82 SE) (4)pOmsk (SAS S IN 7372448): 

Computation of the mean wind distribution 
was made following the standard method em- 
ployed by other authors. At first, heights, 
temperatures, and relative humidities for stan- 
dard isobaric surfaces from 1,000 mb to 200 mb 
were tabulated. Only those soundings that 
reached the 500 mb level were included. Next 
the arithmetic mean of all these elements was 
taken. From the averages of temperature and 
relative humidity, mean virtual temperature 
was computed. Thus a mean sounding for each 
station was obtained and heights were recal- 
culated with the help of Bjerknzs’ tables 
(1910). These computed heights were found 
to be more consistent with hydrostatic equilib- 
rium than the mean of the reported heights. 

Since the flow is fairly zonal above 900 mbs, 
all stations were projected on their own latitude 
to the 76th meridian, and the wind distribution 
was then calculated from the geostrophic wind 
formula. For an isobaric surface this formula 


is U = — (g/f) (>) where U is the zonal 
OY)» 
wind speed taken positive from the west, g the 


acceleration of gravity, f the Coriolis para- 
meter, h the height of an isobaric surface, p the 
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pressure, and y the north-south coordinate 
taken positive toward north. As several of the 
stations are located 1000 to 2000 feet above sea 
level, some difficulty is encountered in the low- 
est levels due to extrapolation of temperature 
below the ground. It turned out to be the best 
to compute the upper wind field using the 
900 mb level as reference surface with zero 
slope, and later add the zonal component at 
ı km as obtained from the mean upper-wind 
charts published by the India Meteorological 
Department. 


Field of zonal motion and absolute vorticity 


The most interesting feature of the wind 
field (fig. 2) is the appearance of three separate 
areas of concentration of west wind. One 
center lies in Siberia, to the north of the Hima- 
layas and the Tibetan Plateau. Its position and 
strength cannot be ascertained completely 
owing to the lack of data above 500 mbs. The 
second center appears immediately to the south 
of the plateau at about 31° N while the third 
lies in the tropics at 15° to 20° N. For the pur- 
pose of future reference these jets will be 
called Siberian, Himalayan, and Equatorial jet 
respectively. 

In order to verify the accuracy of the com- 
puted wind profile a cross section based on the 
upper-wind maps of the India Meteorological 
Department was also prepared (fig. 3). Com- 
parison between the two sections shows a 
more or less identical position of wind maxima 
and minima. However, the wind speeds as 
given by the India Meteorological Department 
are lower, in part presumably owing to the 
short pibal runs on days with strong wind or 
cloudiness. 


Fig. 2. Mean meridional cross-section for January and 
February, 1946 from Bangalore in south to Omsk in 
north along the meridian 76° E. Heavy solid line repre- 
sents the north-south mountain profile. Thin solid lines 
represent constant geostrophic wind (west wind compo- 
nent in mps); dashed lines represent isentropes. 
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Fig. 3. Mean vertical distribution of observed zonal com- 
ponent of west wind (mps) for January along the me- 
ridian 76° E. 


From the computed meridional distribution 
of the zonal wind the field of absolute vorticity 
(fig. 5) has been computed. For currents mov- 
ing approximately along a latitude circle, the 
absolute vorticity (£,) is obtained with little ap- 


proximation from the relation £a = + f— oe 
The most interesting feature of this section is 
the appearance of zero values of absolute vor- 
ticity in a narrow belt 4 degrees latitude wide 
at the southern rim of the Himalayan jet. 
According to the criterion of instability for 
north-south acceleration of air particles moving 
under conservation of absolute angular mo- 
mentum (SOLBERG 1936), little resistance can 
be offered by the stabilizing effect of the 
Coriolis force against southward acceleration 
in this region. Although the validity of this 
criterion for atmospheric motion has not yet 
been demonstrated generally, it will be. of 
interest to look for evidence of descending 
motion near latitude 25° N, in accordance with 
the suggestions of RossBy and COLLABORATORS 
(1947). No explanation for the region of zero 
absolute vorticity, however, is implied since, 
as also pointed out by PALMÉN and NEWTON 
(1948), it is not possible to obtain zero absolute 
vorticity in finite air columns from the opera- 
tion of divergence alone. 


Temperature and potential temperature 


In accordance with the computed wind field, 
a pronounced concentration of the meridional 
temperature gradient (fig. 4) just south of the 
Himalayas is evident. Compared to other 
regions of the world, temperatures encountered 
here in the upper troposphere are remarkably 
low. South of the jet stream region there is an 
actual reversal of the meridional temperature 
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- gradient above 400 mbs between latitudes 25° 
and 18° which, in a mean picture, cannot be 
explained easily on the basis of advective 
processes. 

The configuration of the isentropes (fig. 2) 
conforms closely to the patterns emphasized 
by Rossby and COLLABORATORS (1947). North 
of the jet center we observe a pronounced 
upward bulge of the isentropes and vertical 
crowding in the upper troposphere. On the 
south side, the vertical gradient is small through 
a deep layer of the middle troposphere. This 
stability distribution can be established in a 


more detailed way (fig. 6) by the computation 
90 
az 


where @ is the potential temperature and z 
is the vertical coordinate taken positive up- 
ward. 


of stability (S) from the relation S = O-1 — 


Baroclinity field 


In order to study the distribution of solenoids 
in the mean cross section the average baroclin- 
ity (fig. 7) has been computed. “Baroclinity 
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Fig. 4. Mean temperature and zonal component of geo- 

stropic wind{for Januaryfand February, }1946. Thinfsolid 

lines are soles of westerlyfcomponent?of wind] (mp); 
dashed lines represent isotherms"in 0° 
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Fig. 5. Mean distribution of absolute vorticity (units 

10-5sec-!). Note very low absolute vorticity in a lati- 

tudinal belt of 4 degrees in the high troposphere on the 
southside of the Himalayan jet. 
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Fig. 6. Mean distribution of vertical stability (units 108 
cm"), Note the general decrease of stability downward 
on the southside of the Himalayan jet. 
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is defined as the number of solenoids per unit 
area (N)andis given by — & (=) . Maximum 
ee ey 
baroclinity appears between 400 and 300 mbs 
below and somewhat north of the Himalayan 
jet, but below and somewhat south of the 
Equatorial jet. Two barotropic lines are ob- 
served. One separates Himalayan and Equa- 
torial jets with negative baroclinity above 
300 mbs in correspondence with the decrease 
of the west wind with height between latitudes 
20° and 25°. The other barotropic line lies 
mostly in the transition zone between easter- 
lies and westerlies. 


Discussion of the Himalayan jet 


The new computations bear out the zonal 
wind profile given by Bjerknes in its most 
important aspect—the existence of a jet stream 
center in the high troposphere south of the 
Himalayas. Beyond this point, however, there 
are marked differences, presumably resulting 


fo) 
from the wide variations in longitude of the 
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Fig. 7. Mean distribution of baroclinity (units 108 sec—2). 

Note that the maximum baroclinity is below and some- 

what north of the Himalayan jet while it is below and 

south of the Equatorial jet. Heavy solid lines represent 
zero baroclinity or barotropic lines. 
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stations used by Bjerknes. This observation 
suggests at once that there is no latitudinal 
symmetry around the northern hemisphere. 
A glance at the North American winter section 
of Hess (1948) confirms that this is so, provided 
that the unequal periods used for the computa- 
tion of both cross sections are comparable. In 
order to bring out this asymmetry more 
clearly, the Hess section for 80° W was sub- 
tracted from that for 76° E (fig. 8). 
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Fig. 8. Numerical comparison of zonal component of 

the geostropic wind along 76° E with that along 80° W 

after Hess (1948), the Hess section being subtracted from 

that along 76° E. The anomaly has been expressed in 
mps. 


In the upper troposphere, the flow is mark- 
edly faster near latitude 30° at 76°E com- 
pared to 80° W. Because of more southerly 
position and greater intensity, it is casy to 
explain the difference in the maximum speed 
of the jet stream center. If the southward dis- 
placement of this center at 76° can be ascribed 
to the effect of the Himalayas, its position 
changes little from year to year, and from day 
to day in comparison with regions where a 
permanent obstacle of large size is absent. In 
consequence, higher wind speeds should be 
expected in a mean section at 76° E than at 
80° W where the position of the jet is subject 
to much larger aperiodic variations. 

It is of interest also to compute the mass 
transport in both sections. Total flow of mass 
across 80° W from 1000 to 200 mbs between 
latitude to and 65° N is 570 - 108 tons secu}. 
The latitude belt 15°— 335° contributes 210 : 108 
tons sec! (37 %) to this total. At longitude 
76° E the mass transport is about the same. 
Although it would be necessary to study the 
conditions at many meridians to obtain a con- 
clusive picture, this fact nevertheless suggests 
that mass transport through the troposphere is 
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constant within latitudinal belts of limited 
width. At first glance fig. 8 indicates that 
because of the higher speeds aloft the transport 
should be greater; it turns out, however, that 
there is a compensating effect in the layer 
r000—5$00 mbs. Over India this compensation 
takes the form of the winter monsoonal current. 

It is not possible at this time to pursue this 
reasoning further, although it points to some 
interesting speculative hypotheses concerning 
the dynamical hypothesis of the winter mon- 
soon. It appears permissible merely to state 
that the basic zonal current upon striking the 
plateau splits into two parts that circle the 
mountain ranges. It is suggested that the 
appearance of two west wind maxima is due 
to the effect of the mountain barriers. Some 
qualification, however, is necessary. It is well 
known that on many occasions a marked split 
of the westerly current takes place already in 
Europe, one branch moving along the northern 
rim of the continent and the other travelling 
across the Mediterranean Sea. Although it 
would be difficult to explain the observed 
mean conditions over Central Asia from this 
intermittent feature, it should certainly not be 
denied that on many days there is a jet stream 
center near latitude 30° N extending across 
the entire width of Eurasia. 

Schematically the effect of the Himalayan 
plateau is represented in fig. 9. Here the 
solid line stands for an isobaric surface associ- 
ated with a fast moving current in the initial 
stage before striking the plateau. The dashed 
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Fig. 9. Schematic representation of vertical circulation 

patterns around the plateau. Solid line represents an 

undisturbed isobaric surface associated with a mean jet 

stream center at about latitude 35° N; dashed line repre- 

sents the deformed isobaric surface upon striking the 

plateau. Vertical circulation patterns are indicated by 
the arrows. 
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line represents the deformation of.this isobaric 
surface after the current has passed through a 
portion of the area subject to the mountain 
influence. As demonstrated in the theoretical 
study, by the author (1949), the velocity 
distribution in the currents circling the moun- 
tain ranges will be such as to produce maxima 
in the immediate vicinity of the plateau, 
gradually diminishing on both sides. While* 
this arrangement of current speed can be ob- 
tained in two-dimensional flow up to the 
mountain crest, this approach does not tell us 
why the jet along the southern border con- 
tinues to intensify and sharpen well above the 
mountains. It is necessary to assume that a 
system of vertical circulations exists in the 
upper troposphere. Since hydrostatic equilib- 
rium is always maintained to a high degree 
of approximation, this vertical circulation! 
must be in the sense indicated in fig. 9. 

All supplementary cross sections support 
this contention. The isotherm pattern has 
shown us the concentration of the meridional 
temperature gradient just south of the moun- 
tains and an existence of an area of actual 
reversal of this gradient south of the jet. It is 
especially this latter fact that points to pro- 
nounced subsidence between latitudes 20° and 
25°, since an advective hypothesis is not main- 
tained by the mean section. 

The sections of potential temperature and 
vertical stability bring out the effect that the 
field of horizontal convergence and diver- 
gence necessary to maintain the vertical cir- 
culation cell must have on the stability of air 
column under its influence. It is here assumed 
that this field of divergence has been operative 
for some time before the air reaches the meri- 
dian 76°E, a reasonable contention since the 
plateau of Pamir—the roof of the world— 
lies to the east of this meridian. North of the 
jet center ascent of air can account for the very 
low potential temperature of the high tropo- 
sphere. As the ascent decreases with approach 
to the level of highest wind speed and the 
tropopause, divergence and vertical shrinking 
must occur from continuity requirements. 
Such shrinking is indicated by high values of 
stability in the high troposphere. After the air 
is displaced southward across the jet stream 


1 Such a vertical circulation may be partly due to the 
radiative cooling over the top of the plateau. 
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center, descent and stability decreasing down- 
ward should be encountered. As noted before, 
this stability gradient is clearly observed. This 
description accords with that given by RossBy 
and COLLABORATORS (1947). 

The observed gradient of absolute vorticity 
shows that the vertical motion pattern as out- 
lined is at least dynamically possible, even 
though not proved, in view of the assumptions 
on which the Solberg instability criterion is 
based. However the baroclinity section offers 
strong support in favor of indirect circulation. 
If the jet stream were maintained through 
direct circulation, we should encounter maxi- 
mum baroclinity south rather than north of 
the jet stream. 

As final evidence we may turn to rainfall 
conditions during winter. If the Himalayan 
jet is orographically produced, its position 
should suffer little change from year to year 
as stated before. In consequence a zone of 
minimum rainfall should occupy central India 
during winter. This is born out by actual con- 
ditions (fig. 10). Noteworthy also is the narrow 
rainy zone extending along the northern edge 
ot the sub-continent where “winter” de- 
pressions move eastward underneath the jet 
stream. 

It is not possible at this time to propose a 
concrete explanation of the vertical circulation 
system. It may, however, be suggested that 
the surface frictional stresses of the enormous 
plateau area can be an important factor. This 
suggestion is backed by the fact that the moun- 
tains are extremely rough, and that they rise 
steeply along their western border from low 
elevations to the vicinity of the soo mb surface 
or higher, so that the upper-air current must 
strike them at a high speed. 


Discussion of the Equatorial jet 


The structure of atmosphere over southern 
India is very different from that over the 
northern part of the country. Here we find 
a region of west wind concentration over a 
lower layer of casterlies. Inspection of tem- 
perature, potential temperature and baroclinity 
fields (figs. 2, 6 and 7) indicates that there are 
enough solenoids to maintain both lower 
easterlies and upper westerlies, provided that 


.there is a vertical circulation system in the 
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work-producing sense. Although the existence 
of such a circulation has been assumed by 
writers on the general circulation for many 
years, it was not until quite recently that the 
actual strength of inflow into the equatorial 
zone and the distribution of convergence in 
that region has been computed (1949). It 
follows from these computations, together 
with the continuity principle, that the previ- 
ously assumed direct solenoidal circulation in 
low latitudes should be a reality. 

This argument also finds support in the 
baroclinity field. In pronounced contrast to 
conditions over northern India maximum 
solenoidal concentration in the equatorial zone 
appears below and south of the equatorial jet. 
Direct vertical circulation should build up a 


region of west wind concentration north of 


the maximum as is observed. 

In consequence, ascending motion should 
occur south of the jet stream center and 
descent to its north. Winter rains should be 
confined largely to the region south of the jet, 
as is indeed the case (fig. 10). It therefore 
appears permissible to conclude that the 
equatorial jet over Indo-Pakistan can be con- 
sidered as an example of the “HADLEY” type 
(1735) of circulation. 


Fig. 10. Mean distribution of rainfall over Indo-Pakistan 
in January in inches. 
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Abstract 


The percentage variations of twenty-four hour means of C. R. in Stockholm (magn. lat 58°) 
amount to + 5% of the mean value over a quiet period of two months’ duration. The two 
atmospheric effects, the mass absorption due to ionisation and the meson decay, contri- 
bute + 4%. The residual variations, in excess of statistical fluctuations of 0.1%, depend 
upon other factors, such as geomagnetic activity and solar influences. 

The true absorption coefficient for zenith radiation (f) is calculated to have a value of 
1.27: 107$ cm?/g and the rate of decay (P’) is found to be about 0.060 pr km. The meson 
producing layer is attributed to the 200 mb surface, the highest level reached by the daily 
radiosonde ascents in Stockholm. 

The most reliable decay value is obtained by considerations of the thermal structure of 
air-mass changes in travelling deep depressions. 


Introduction 


I am indebted to Dr Malmfors for copies of the 
continued counter records from Stockholm (MALMm- 
FORS 1949) which have been used in the investiga- 
tions discussed in this paper. The first objective, of 
an analysis of the intensity values Z, N and S of 
C. R., was to deduce the true absorption and decay 
coefficients of particles from the three central direc- 
tions with the mean zenith distances 9 = 0°, 35° N 
and 35°S. 

The change dN in the number of particles pe- 
netrating from a direction deviating 9° from the 
zenith depends upon ionisation and decay and can 


be expressed by 
dH = — NßdB sec 9 — NB'dH sec à 


where dB in mm Hg is the change in the air-mass 
over the station and dH the height-changes in m 
of the meson producing layer in the lower stra- 
tosphere ß representing the true absorption coefli- 


I 
cient and ß’ the mean rate of decay. p’ = = L the 


mean range of the meson before disintegration. 


With the mean values N, B and H for the whole 
period in question we have 


N—N =— NB (B—B) sec 9 — NB (H—H’) secd 


The daily mean-values of the airmass, the baro- 
metric pressure, were taken from the hourly values 
of the Meteorological station at Stockholm- 
Bromma airport. 

Tre results of counter-soundings to upper levels 
(PFOTZER 1936) show the meson producing layer 
to be at a pressure of about 100 mb. Unfortunately 
the daily aerological and wind soundings at Stock- 
holm generally do not reach higher than the 200 
millibar surface. These soundings, which have been 
employed in the present investigation, are made 
twice daily at 5° and 15h GMT. The geopotentials 
of the 200 millibar surface are taken from these 
soundings and the corresponding height values have 
been given in geometric metres. 

The analysis of the daily means has given the 
following regression coefficients y and y’ in the 
formula 
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Fig. 1. Observed percentage variations of the daily 

means of Z, N och S intensities, the barometrique 

pressure at the ground and the height of the 200 mb 
surface. 


N—N =—u (B—B)—w’ (H— HM): 


Z — 25 423 = — 44.05(B — 757.3) — 1.629(H — 
TI 340) 

N — 16 263 = — 29.94(B — 757.3) — 0.920(H— 
II 340) 

S — 15 645 = — 28.45(B — 757.3) — 0.998(H— 
II 340) 


As shown in fig. 1 the observed variations of the 
daily means of Z, N and S amount to about + 5% 
of the means for the whole period. However they 
are reduced to + 1 % by considering both mass 
absorption and decay (the thick curves in fig. 2). 
The thin curves in fig. 2 represent the percentage 
deviation after reduction to a constant pressure 
of 757.3 mm Hg. Other factors are responsible for 
residual changes larger than the statistical Auctua- 
tions of about 0.1 % of the mean intensity. Though 
not yet fully elucidated the influence of magnetic 
activity and of solar influences on the primary 
radiation seem to be of particular importance. 


1. Barometric coefficient. We obtain the true 
absorption coefficient ß from the barometric 
regression coefficient à by means of the equation 


b= NB sec Ÿ 
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With the above values of # the absorption coeffi- 
cients Boe and Bas for zenithal and for an incidence 
of 35° are 


Boo = 1.73 % per cm Hg or 1.27x 1073 cm?/g. 


Paz: = 


This result is in accordance with our knowledge of 
the C. R., for with a greater angle of incidence the 
weaker component of the C.R. is absorbed by 
the greater air-mass in the oblique direction. 

As at sea level the hard meson-component is 
about 75 % of the total zenithal radiation and as 
the N or S intensities are found to be in about the 
same relation to the Z intensity, we may suppose 
that the oblique radiation N or S is about equivalent 
with the hard meson-component of the zenithal 
radiation. The absorption coefficient of this meson- 
component must then be of the order of 1.10 x 1073 
cm?/g. 

The calculated values of B are lower than the 
BE-value LINDHOLM (1944) deduced without ac- 
counting for the disintegrating process, as H is also 
a function of B according to the hypsometrical 
formula. 


150.5 » » DO TRIO TOM cm?/g 


2. Decay coefficient. From the regression coefficient 
w we deduce the rate of decay 


for. Z Bo = 6.5% per km and 
» N+S B'as — 49% » km. 


Considerations of the thermal structure of travel- 
ling cyclones, causing different slopes of the isobaric 
surfaces at upper and lower levels, suggest a method 
for determination of the rate of decay independently 
of the barometric effect. During the period in- 
vestigated there were two deep depressions passing 
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Fig. 2. Percentage variations of daily means of Z, N 
and S. Thin curves represent variations after reduction 
to the mean pressure (757.3 mm Hg). Thick curves 
represent variations after reduction both to mean 
pressure and mean height of the 200 mb. surface. 
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over Scandinavia, one during the days s—7 January 
1948 with a pressure minimum at Stockholm of 
730.6 mm Hg at 16 on the 6th and another one 
during the days 12—14 January with a minimum 
of 735.5 mm Hg between o—3h on the r4th. 

Fig. 3 shows the average pressure variation at the 
ground and furthermore the variation of the height 
of the 200 mb surface and of the Z and N+ S 
values during the passages of the two depressions. 
The pressure and height values are taken at three 
hour intervals before and after the occurance of the 
pressure minimum at the ground. The values of 
Z and N + Sare the sum of the number of coincid- 
enses at the same time interval on the two occasions. 
The points indicate the values of Z and N+ S$ 
reduced to a mean pressure of 737,5 mm Hg. 

The corresponding observed mean values of H 
(200 mb) in meters at a mean pressure of 737.5 mm 
Hg before and after the passage of the depression 
and the sum of Z and N + S at the two events are 
given in Table 1. 


Table 1. Mean values of Z, N + S and height 
of the 200 mb surface for constant pressure at 
the ground during depression passages. 


Hour 
from 
S 
barometer mb 
Sy t the : 
minimum surface 
passage 


66,460 
66,820j— 2.93 


12%, before|/737°5|11,247 
1214 after |737.5|11,124 


53,000 
53,400|— 3.25 


From the value of w the rates of decay are calculated 


B'; = 0.0061% pr meter or 6.1% pr km 


B’nis = 9.0036 % » » DE | 


If a great number of cyclone passages are treated in 
this manner the determination of f’ can be made 
with great accuracy without knowing the baro- 
metric effect. 

In fig. 4 the values of Z and N+ S reduced to 
737.5 are represented as functions of Hand H- sec 35° 
resp. The correlation coefficient between Z and 
N + S reduced to constant pressure and H?% mb 
is g = — 0.92. By means of the least square 
method we obtain: 


Fig. 3. Average variations of B, H, Z and H+S 

during passages of two deep depression. (The points 

indicate values of Z and N + S reduced to a mean 
pressure of 737.5 mm Hg.) 
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The high correlation between the intensity 
reduced to constant pressure and the height of the 
200 mb surface support the assumption that the 
variations in height of the 200 mb surface differ 
only insignificantly from those of the meson 
producing layer assumed to be at about 100 mb. 

SCHERHAG’S (1948) northern hemisphere charts 
of the mean absolute topography of the upper 
troposphere and the lower stratosphere, drawn by 
means of soundings 1941—44, show but little change 
from the 225 to the 96 mb surface of the difference 
between the extreme values equator (max) — pole. 
Expressed in geodynamic decameters (gdm) this 
variation is found to have the following values: 


Level | Equator Pole | Difference 
ape 1] 225 mb |1,136 gdm| 980 gdm| 156 gdm 
“HT 906 9 27,04326) 217,492 0) 156 » 

G4 Sy 7.008 > 6 » 
I { ER wa) al it ; 9 
July 96 » |1,688 » [1,624 » 64 » 


It is interesting to note that, if we take the corres- 
ponding height values for the 300 mb level, we 
find a very little change in the calculated value 
namely B’ = 5.92% per km. 

KOLHÔRSTER and MATTHEW (1930) found ff, = 
= 113X10 3 cm?/g by recording double coin- 
cidences at Potsdam, and compensating for the 
variation of atmospheric pressure with a variable 
wood screen. This value corresponds to f’y = 6.6 ° 
per km and is quite close to that given above. 

From a two months record in London in 1943 
DUPERIER (1944) found, from threefold coincidences, 


a lower value Bo = 5.4 % per km at a mean pressure 
of 756.7 mm. This value was deduced by means of 
height values of the 100 mb surface over southern 
England. The LE na cocfficent, however, was 
id to be 2.28 9 o pr cm Hg so that the sum of 
the two en is the same in both cases. 

It is interesting to find how well the height of the 
200 mb surface can be determined by records of 
the cosmic radiation. With the decay formula 
given above the height values have been calculated 
from the intensity values of Z reduced for pressure 


e (H) = + \ Se 


the individual height values is found to be + 32 mi.e. 

3 pro mille or not tar from the accuracy of x actual 
soundings. This result was obtained by averages from 
the two passing cyclones. 

For the whole material we find e (H) = 
or 8 pro mille. 

3. Temperature coefficient. The temperature coeff- 
cient is obtained from 


variations. The mean error 


+ 96 m 


Ber. dH 
OBER NES de 


If we define the mean temperatur T,, of the air 


m 
column between B and the 200 mb surface by the 
hypsometric formula 


Bo 
Ty log — km 
>? 


H500 mb = 0.0674 - 


> IT a Le 
I OIO mb = 0 0474 and 


with Bp’ = 0.06 per km we have & = —0.06X 
X 0.0474 = — 0.0028 pr °C or & = — 0.28% pr °C. 


we find, for Bo = 


We obtain for the air column between B and 


dH £ 
100 mb ( ) = 0.00677 and assuming the same 


m 


value fp’ = 0.06 pr km 
& = — 0.00406 or — 0.41% per °C. 


DUPERIER (1944) has found the value of « = —0.35 % 
per °C. Other methods for calculation of the tem- 
perature effect have given values of « between 
— 0.30 and — 0.40 pr °C. 


4. Meson-decay and latitude effect of cosmic radiation. 
The mean slope equator — pole of the isobaric 
surfaces in the upper troposphere and lower strato- 
sphere, and also of the layer where mesons are 


- 
u 
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generated, must have some bearing on the inter- 
pretation of the latitude effect of C. R., which is 
found by Cray (1932), HoERLIN (1933) and in the 
extensive observations of A. H. Compron (1933). 
From the fact that a latitude effect — an increase 
of the intensity of 12 per cent of the radiation —, 
is found, when going from the equatorial region 
to 55 °N geomag. latitude, it is concluded by these 
authors that this part of the radiation consists of 
corpuscles deviated by the earth’s magnetic field 
before entering the earth’s atmosphere. 

In table 2 the meridional changes of the topo- 
graphy of the 96 mb surface are given in dyn deca- 
meter (gdm) as taken from SCHERHAGS mean 
chart for January and July. The great seasonal 
variation of the slope of the 96 mb surface with 
maximum in the wintermonths and minimum in 
summer is evident from this table. 


The decay causes a variation with latitude of the 
mean intensity of C. R. as is shown in column 3 and 
5 of table 2 for January and July respectively. 

It is thus evident that in the winter months the 


Table 2. Average absolute topography of the 
96 mb surface over the northern hemisphere ac- 
cording to SCHERHAG (1948) 


January July 

ats no 

Mean Mean Rel. Mean | Mean Rel. 

Geop. gdm.!Int. C.R. Z Geop. gdm.|Int. C.R. Z 

(SSeS emp Kr 370 A 
90° N 1,496 1.097 1,632 1.0477 
80 1,504 1.088 1,636 I.O15§ 
70 1,540 1.058 1,640 1.012 
60 1,564 1.049 1,644 1.010 
so 1,580 1.040 1,652 1.00$ 
40 1,596 1.033 1,670 0.994 
30 1,610 1.027 1,682 0.987 
20 1,624 TOR 1,676 0.990 
10 1,640 1.005 1,664 0.998 
fe) ' 1,648 1.00 1,660 1.00 


decay effect contributes with about one third in 
the latitude effect, up to 55° latitude. Still farther 
north the meson decay must cause further increase 
of the average intensity amounting to 10 % in the 
polar region. 
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BOOK REVIEW 


B. GUTENBERG AND C. F. RICHTER, Seismicity 
of the Farth and Associated Phenomena. Prin- 
ceton Univ. Press, Princeton, New Jersey, 1949, 


273 PP: 


When two of the world’s most distinguished 
scismologists have been working on a subject for 
more than ten years, one may expect important 
results to come out. These expectations are more 
than satisfied by this book, which contains the 
important results obtained at the Seismological 
Laboratory at Pasadena. The book is epoch-mak- 
ing in geophysics, particularly because it treats 
a statistically homogeneous material. This is facil- 
_itated by a consequent use of Richter’s magnitude 
scale for earthquakes, originally introduced in 
1935, and considerably developed in subsequent 
years by Gutenberg and Richter. Its importance 
lies in the fact that the magnitude is strictly defin- 
ed in instrumental quantities, thus completely 
objective. 

The reviewer is very much impressed by the 
large amount of work which lies behind and not 
less by the accuracy and conscientiousness with 
which it has been carried out, e. g. all the deter- 
minations taken from the International Seismol- 
ogical Summary have been carefully examined, as 
well as all other previous determinations. All 
seismologists in the world must be very grateful 
to the authors that their determinations have been 
compiled in the tables in the book, pp. 117— 
252, containing date, time, location, depth, magni- 
tude etc. for about 3000 earthquakes, a list of the 
greatest value for seismological investigations. Pp. 
253— 267 contain a table of about 450 active 
volcanoes. Another feature of great help in sub- 
sequent investigations is the large list of references 
(pp. 104—116), containing about 560 numbers, 
including not only papers on seismological studies 
of various earthquakes (both instrumental and 
directly in the carthquake regions), but also on 
vulcanology, gravity, geology, topography etc. 

After a few introducing pages we find an 
interesting chapter on frequency and energy of 
earthquakes (pp. 16— 25). It gives statistical results, 
which thanks to the magnitude scale are more 
reliable than hitherto. The basic relations are one 
for the energy E corresponding to the magnitude 
M and one for the annual number N of earth- 
quakes of magnitude M: 


10Log E = 12 + 1.8 M 
Log N=a+b (8 —M) 


a, b are constants, dependent on the depth of the 
focus. We learn that the total number of true 
carthquakes may well be of the order of a million 
each year; that the mean annual energy released 
in earthquakes is about 1027 ergs, and that the 
few large earthquakes contribute most to the 
energy. Table 7 on p. 22 gives statistics con- 
cerning the geographical distribution of earth- 
quakes. Among other things we find that of the 
shallow earthquakes (depth = 60 km) about 80 % 
belong to the Circum-Pacific Belt, of intermediate 
earthquakes (depth 70—300 km) 91.3% belong 
to the Circum-Pacific Belt and the rest to the 
Transasiatic Belt, and that all deep earthquakes 
(depth > 300 km) are located to the Circum- 
Pacific Belt. 

The greatest part of the text (pp. 28—94) is 
devoted to the regional discussion, where region 
by region round the earth is treated with regard 
to earthquakes, volcanoes, gravity anomalies, geo- 
logy etc. The Pacific arc structures have been 
studied with particular care. A vertical cross-sec- 
tion of a typical Pacific arc is given in Fig. 6, p. 
29. The book contains almost exclusively facts 
and very few hypotheses. But in spite of the 
abundance of facts it is written in an enjoyable 
form. The accompanying maps are excellent; one 
could possibly wish that they would contain 
more of the geographical names given in the text. 
The facts and the results are too many to be 
mentioned here. One chapter (pp. 94—97) treats 
tsunamis (seismic sea waves) and contains many 
references. Pp. 97—101 are devoted to »Mechan- 
ism», and contain several interesting ideas. A ques- 
tion of interest is naturally to what extent recor- 
ded earthquakes for somewhat more than 40 years 
are representative, in other words it would be 
very interesting to compare the earthquake distri- 
bution given here with that of another 40—so 
years. 

I hope to have given by this review the im- 
pression that the book is of the greatest value for 
all who have something with the earth to do in 
their science, geographers, geologists, geophysicists 
of all kinds, particularly seismologists. It is hoped 
that the book will be a standard work in its field 
for many years, but it is also hoped that it will 
be continued as time goes on and more data are 
accumulated. 


Markus Bath 
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